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Abstract 


The  intense  fluctuating  pressures  observed  in  open  cavities  are  extremely  detrimental  hence  the 
ability  to  control  the  flow  to  minimize  them  is  an  important  challenge.  Many  open  loop 
methodologies  have  proven  effective  at  this  for  single  operating  conditions  however  they  often 
are  detrimental  at  off  design  conditions.  Closed  loop  flow  control  is  applicable  over  a  range  of 
flow  conditions  and  is  the  thrust  of  many  current  and  past  studies.  However  in  order  to  advance 
its  effectiveness  there  is  a  need  for  a  better  understanding  of  the  flow  physics  and  the  effect  that 
control  has  on  it.  This  study  represents  a  step  in  the  direction  of  developing  the  tools  that  will 
allow  for  a  better  understanding  of  this  flow  physics  along  with  an  improvement  adaptive  control 
methodologies.  One  aspect  presented  in  this  study  involves  the  application  of  an  estimation 
methodology  to  obtain  the  time  dependent  flow  properties  in  high  speed  applications  and  then 
calculate  the  relevant  flow  properties  associated  with  generation  of  the  pressure  fluctuations.  In 
conjunction  with  gaining  insight  into  the  flow  physics  we  have  also  developed  improved 
actuators  and  an  adaptive  flow  control  methodologies  which  demonstrate  reductions  in  the  broad 
band  and  tonal  components  of  the  fluctuating  surface  pressures.  In  addition  wc  present  a  study 
which  was  conducted  using  successful  open  loop  control  to  demonstrate  our  ability  to  determine 
the  effects  of  the  control  on  the  underlying  flow  physics. 
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1 


Introduction 


The  primary  goal  of  this  research  is  to  develop  a  better  fundamental  understanding  of  the 
sources  generating  the  large  surface  pressure  fluctuations  in  cavity  flows  so  that  more  efficient 
physics  based  control  can  be  used  in  the  future.  This  first  chapter  introduces  the  flow  physics 
of  cavity  flows  by  way  of  a  brief  discussion  then  we  will  discuss  some  of  the  previous 
applications  of  control  to  cavity  flows.  This  is  followed  by  a  brief  statement  of  the  objectives  of 
the  study  which  is  in  turn  followed  by  section  discussing  the  layout  of  the  rest  of  this  document. 

1.1  Motivation 


From  the  limited  review  that  will  be  presented  below,  it  is  apparent  that  while  progress  has 
been  made  towards  effective  active  suppression  concepts,  we  do  not  yet  understand  the  key 
physical  meehanisms  that  need  to  be  leveraged  in  the  design  of  an  effective  control  system 
without  iteration.  The  current  effort  is  our  initial  effort  to  develop  adaptive  control  concepts  and 
the  tools  which  will  allow  us  to  investigate  their  effects  on  the  flow  physics.  Through  this  and 
follow  on  efforts  we  hope  to  provide  a  methodology  that  will  us  to  develop  better  control 
methods  along  with  having  a  better  understanding  of  the  flow  physics  effected  by  the  control. 


1.2  Background 

In  this  section  we  will  review  the  basic  physics  of  flow  over  an  open  cavity  along  with  a 
detailed  review  of  many  studies  which  have  used  flow  control  (open  and  elosed  loop)  in  an 
attempt  to  reduce  fluctuating  surface  pressure.  This  review  is  important  as  it  provides  the  back 
drop  of  what  has  been  previously  accomplished  and  set  the  stage  for  the  work  presented  here  and 
follow  on  efforts  to  understand  the  effects  of  control  on  the  flow. 


1.2.1  Cavity  Flow  Physics 

Flow  over  an  open  cavity  has  important  characteristics  that  are  representative  of  those  in 
many  industrial  configurations  as  well  as  exhibiting  important  fundamental  phenomena.  Under 
eertain  geometrical  configurations  the  flow'  is  known  to  be  dominated  by  a  large  amplitude  self- 
sustained  oscillations  that  originate  from  a  complex  aeroacoustic  feedback  process.  This  process 
is  dominated  by  the  growth  and  eonveetion  of  instabilities  in  the  shear  layer  at  the  cavities 
opening  and  how  this  shear  layer  interacts  with  the  cavities  rear  wall.  The  oscillations  of  the 
shear  layer  and  the  impact  of  the  free  stream  flow  on  the  aft  wall  of  the  cavity  serve  as  a  source 
which  fulfill  the  self-sustained  oscillations  and  result  in  large  pressure  fluctuations  throughout 
the  cavity.  These  large  pressure  fluctuations  are  undesirable  in  most  applications  and  the  need 
for  a  better  understanding  of  this  process  so  that  a  means  to  alter  it  can  be  developed  is  still  a 
topic  of  many  ongoing  studies  even  after  over  a  half  century  of  research  on  resonating  cavity 
flows. 
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Flow  over  open  cavities  has  been  studied  over  the  past  several  decades  from  several 
different  perspectives.  The  original  studies  which  slanted  towards  aircraft  applications  arc 
generally  attributed  to  Roshko  (1955)  and  Krishnamutry  (1955).  The  more  than  50  years  of 
studies  of  flows  over  open  cavities  which  have  followed  have  led  to  many  advances  in  the 
understanding  of  the  important  flow  features  of  cavity  flows  and  the  governing  parameters  which 
dictate  its  behavior.  The  early  studies  revealed  an  unsteady  momentum  flux  into  and  out  of  the 
cavity,  which  results  in  a  "breathing"  or  "pumping"  action.  This  pumping  action  causes  the 
stagnation  streamline  of  the  shear  layer  to  fluctuate  up  and  down  the  aft  wall  of  the  cavity,  which 
has  been  widely  thought  to  be  responsible  for  sound  production  in  the  cavity.  The  study  of 
Rossitcr  (1964)  is  regarded  as  one  of  the  significant  benchmark  cavity  flow  studies  and  resulted 
in  a  semi-empirical  relationship  that,  with  slight  modification,  predicts  the  frequencies  of  the 
resonant  modes  in  open  cavities.  Rockwell  and  Naudascher  (1978)  provided  an  excellent  review 
of  the  work  conducted  on  cavity  flows  through  the  late  1 970's.  One  of  the  more  comprehensive 
studies  was  that  of 

Heller  and  Bliss  (1975),  which  was  an  analytical  and  experimental  program  that 
examined  cavities  exposed  to  both  subsonic  and  supersonic  free  stream  flows.  In  their  study  they 
postulated  mechanisms  of  how  the  shear  layer  and  free  stream  flow  interact  with  the  aft  wall  of 
the  cavity  and  also  evaluated  several  concepts  for  suppressing  the  discrete  tones,  such  as  slanting 
the  aft  wall.  A  water  table  was  used  to  visualize  the  aforementioned  interactions,  and  the  authors 
claimed  that  the  important  mechanism  is  not  necessarily  the  interaction  of  the  shear  layer  with 
the  aft  wall  but  instead  the  free  stream  flow  impinging  on  the  aft  wall  when  the  shear  layer  enters 
the  cavity.  More  recently,  schlicrcn  images  of  a  cavity  in  transonic  flow  by 

Heller  and  Delfs  (1996)  have  largely  substantiated  the  earlier  proposed  mechanisms.  In 
addition,  there  have  been  many  studies  that  have  concentrated  on  collecting  databases  for 
characterizing  parametric  effects,  such  as  Tracy  and  Plentovich  (1997)  and  Dix  and  Bauer 
(2000).  From  a  review  of  the  previously  mentioned  studies  and  many  others,  too  numerous  to 
list  here,  one  can  make  several  observations.  The  first  is  that  the  features  of  cavity  flow  are 
largely  determined  by  the  coupling  of  the  shear-layer  instabilities  with  the  acoustic  properties  of 
the  cavity  and  it  is  the  interplay  between  the  shear-layer  instabilities  and  the  cavity  acoustics  that 
cause  each  to  excite,  modify,  and  amplify  the  characteristics  of  the  other. 

Secondly,  the  impingement  of  high  speed  flow  on  the  aft  cavity  wall  results  in  a  pressure 
fluctuation  on  the  cavity  boundary  that  causes  the  downstream  reattachmcnt  region  to  act  as  an 
acoustic  source. 

This  pressure  source  is  part  of  a  feedback  loop  between  itself  and  the  shear-layer  via  its 
interaction  with  the  the  shear-layer.  While  shear-layer  instabilities  and  wall  pressure  fluctuations 
exist  for  all  flow  conditions  over  a  cavity  of  sufficient  size,  their  existence  docs  not  guarantee 
this  resonant  condition. 

The  existence  of  this  resonant  condition  is  governed  by  many  quantities,  such  as  the 
length-to-depth  ratio  (L/D),  cavity  volume,  ratio  of  boundary  layer  thickness  to  cavity  length 
(5/L),  Mach  number  of  the  free  stream  flow,  ratio  of  shear-layer  convection  velocity  to  free 
stream  velocity  (k),  temperature  ratio  between  the  cavity  and  the  free  stream,  and  the  state  of  the 
incoming  boundary  layer. 

With  the  advent  of  more  advanced  tools  for  studying  fluid  dynamics  (P1V  and  time 
dependent  simulations)  there  have  been  several  recent  studies  which  have  looked  at  specifics  of 
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the  flow  field  for  a  wide  variety  conditions.  In  the  existing  literature  studies  which  involve 
quantitative  flow  measurement  appear  to  be  for  cavities  in  incompressible  flows  (e.g.  Grace  ct  al, 
2004,  Ashcroft  &  Zhang,  2005,  Ukeiley  et  al,  2005  and  Kang  et  al,  2008)  or  supersonic  flows 
(such  as  Zhuang,  et  al,  2006  and  Unalamis  et  al,  2001)  leaving  the  gap  for  the  data  in  this 
manuscript  to  fill  in.  These  studies  have  reinforced  the  importance  of  the  incoming  boundary 
layer  and  have  allowed  for  studying  the  vortical  structures.  For  example  Grace  ct  al.  (2004) 
studied  the  effects  of  turbulent  and  laminar  boundary  layers  on  the  flow  features  of  an  length  to 
depth  ratio  4  cavity  in  an  incompressible  flow.  They  found  many  of  the  flow  properties  to  be 
altered  by  the  state  of  the  boundary  layer  and  asserted  that  the  in  order  to  predict  whether  or  not 
resonance  might  occur  or  match  numerical  data  it  is  crucial  to  match  the  conditions  of  the 
upstream  boundary  layer.  The  lower  speed  such  as  Kang  et  al.  (2008)  and  Aschcroft  and  Zhang 
(2005)  along  with  the  limited  transonic  studies  such  as  Hirahara  et  al.  (2007)  which  have 
acquired  velocity  measurements  have  been  able  to  clarify  some  of  the  properties  of  the  shear 
layer  and  its  interaction  with  aft-wall  as  was  originally  proposed  in  Rockwell  and  Knisley 
(1979).  Additionally  these  studies  have  show  ed  linear  growth  of  the  shear  layer  width  and  hence 
the  turbulent  length  scales  in  the  shear  layer  along  with  some  descriptions  of  the  randomly 
sampled  instantaneous  velocity  and  vorticity  fields  demonstrating  the  large  scale  structures  in  the 
flow.  Murray  and  Ukeiley  (2003)  developed  a  novel  technique  for  estimating  the  time 
dependence  of  the  large  scale  structure  and  applied  it  to  the  data  being  presented  in  this  work  in 
Murray  and  Ukeiley  (2007).  There  it  was  shown  how  the  propagation  of  the  large  scale 
structures  in  the  shear  layer  is  accompanied  by  a  minima  in  the  surface  pressure  and  the  idea  of 
the  number  of  turbulent  structures  across  the  cavities  opening  being  associated  with  which 
Rossiter  mode  is  dominant  was  validated. 


1.2.2  Control  of  Cavity  Flows 


Techniques  to  suppress  cavity  oscillations  can  be  classified  in  several  ways.  In  this 
paper,  we  choose  the  classification  detailed  in  of  Cattafesta  et  al  (2008)  to  be  consistent  with 
terminology  used  in  active  noise  and  vibration  control.  Active  control  provides  external  energy 
(e.g.,  mechanical  or  electrical)  input  to  an  adjustable  actuator  to  control  the  flow,  while  passive 
control  techniques  do  not.  Passive  control  of  cavity  oscillations  has  been  successfully 
implemented  via  geometric  modifications  using,  for  example,  fixed  fences,  spoilers,  ramps, 
(Heller  &  Bliss,  1975  and  Shaw,  1979)  and  a  passive  bleed  system  (Chokani  and  Kim,  1991). 
Note  that  some  control  devices  considered  passive  by  this  classification  extract  energy  from  the 
flow  itself  and  have  been  called  ‘ active ’  by  other  researchers.  Pertinent  examples  include 
unpowered  or  passive  resonance  tubes  (Stanek  ct  al,  2000)  and  cylinders  or  rods  placed  in  the 
boundary  layer  near  the  leading  edge  of  the  cavity  (Mcgrath  &  Shaw,  1 996  and  Ukeiley  et  al, 
2004).  These  devices,  described  further  in  Section  3,  are  sometimes  referred  to  as  active  because 
they  provide  an  oscillatory  input  to  the  flow,  but  their  effect  on  the  flow  cannot  be  adjusted 
without  either  changing  the  flow  conditions  or  changing  the  device  itself. 

Active  control  is  further  divided  into  open-  and  closed-loop  approaches.  By  its  very 
definition,  closed-loop  control  implies  a  feedback  loop,  in  which  some  flow  quantity  is  directly 
sensed  or  estimated  and  fed  back  to  modify  the  control  signal  (Distcfano  et  al,  1990).  Open  loop 
corresponds  to  the  case  when  there  is  no  feedback  loop. 
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A  further  non-standard  but  useful  classification  of  closed-loop  flow  control  is  that  of 
quasi-static  vs.  dynamic  feedback  control.  The  quasi-static  case  corresponds  to  slow  tuning  of  an 
open-loop  control  approach  and  occurs  when  the  time  scales  of  feedback  arc  large  compared  to 
the  time  scales  of  the  plant  (i.e.,  flow).  This  approach  is  particularly  relevant  in  nonlinear  fluid 
dynamic  systems,  where  the  fundamental  notion  of  frequency  preservation  in  a  linear  system 
does  not  hold.  The  quasi-static  approach  was  successfully  used  by  Shaw  and  Northcraft  (1998). 
The  usual  dynamic  compensation  case  corresponds  to  the  situation  when  the  above  time  scales 
are  commensurate.  This  can  be  implemented  using  an  analog  (see,  for  example,  Williams  ct  al, 
2000)  or  “real-time”  digital  control  systems  (Cattafcsta  ct  al,  1999).  In  this  context,  “real  time” 
refers  to  the  situation  in  which  the  control  signal  is  updated  at  the  sampling  rate  of  the  data 
system,  and  the  actuator  responds  to  the  flow  state  changes  at  the  time  scales  of  the  dynamics. 


1.2.2. 1  Open  Loop  Suppression  Techniques 

It  is  impossible  to  include  all  prior  research;  there  arc  numerous  other  (mostly  passive) 
studies  that  have  not  been  included  here  due  to  size  constraints.  Some  key  observations  are 
discussed  here. 

Sarohia  and  Massier  (1977)  studied  the  efficacy  of  steady  mass  injection  at  the  base  or 
floor  of  two  different  axisymmctric  cavity  models  for  both  laminar  and  turbulent  boundary 
layers.  While  base  injection  was  effective  at  suppressing  cavity  tones,  large  mass  flow  rates 
were  required  Bc  of  5-15%.  The  value  of  that  will  be  used  through  out  this  section  involves 

normalizing  the  mass  flow  rate  of  the  actuator  normalized  by  the  product  of  free  stream  density, 
free  stream  velocity  and  the  area  of  the  cavity. 

Samo  and  Franke  (1994)  studied  static  and  oscillating  fences,  and  also  steady  and  pulsed 
injection  (at  0°  or  45°  with  respect  to  the  free  stream  direction)  at  transonic  speeds  near  the  cavity 
leading  edge.  Blowing  coefficients  B  of  up  to  7%  were  used.  While  the  static  fences  provided 

the  best  suppression,  the  bandwidth  of  the  mechanical  fences  was  limited  to  <  220  Hz,  while  the 
pulsed  injection  was  <  80  Hz.  These  frequencies  were  at  least  an  order  of  magnitude  lower  than 
the  frequencies  of  the  cavity  tones,  and  thereby  constituted  a  quasi-static  or  low-frequency 
forcing.  Nonetheless,  they  represented  an  important  step  in  the  evolution  of  active  control  of 
cavity  oscillations,  both  in  terms  of  approach  and  also  the  introduction  of  scaling  laws  for  such 
actuators. 

Vakili  and  Gauthier  (1994)  obtained  significant  acoustic  tone  attenuation  with  steady 
normal  mass  injection  through  variable-density  porous  plates  upstream  of  the  cavity  leading  edge 
at  Mach  1.8  using  Bc  ~  4%.  They  attributed  the  attenuation  to  a  thickening  of  the  cavity  shear 
layer  and  a  corresponding  alteration  of  its  instability  characteristics. 

McGrath  and  Shaw  (1996)  subsequently  studied  mechanical  oscillations  of  hinged  flaps 
at  frequencies  up  to  35  Hz  over  a  range  of  subsonic  and  supersonic  Mach  numbers.  Similar  to 
the  Samo  and  Franke  experiments,  the  forcing  frequency  was  an  order  of  magnitude  lower  than 
the  resonant  tone  frequency.  The  static  and  oscillatory  deflections  were  on  the  order  of  the 
boundary  layer  thickness  8  and  were  shown  to  be  effective  despite  their  limited  bandwidth. 


7 


MeGrath  and  Shaw  were  the  first  to  study  the  effeet  of  a  cylinder  plaeed  in  the  upstream 
subsonie  boundary  layer.  Because  of  the  well-known  shedding  eharaeteristies  of  a  eireular 
cylinder  of  diameter  d ,  St  =  fd/U  *  0.2,  over  a  wide  Reynolds  number  range,  this  deviee  was 
called  a  high-frequeney  tone  generator  (HFTG).  The  cylinder  was  capable  of  producing 
substantial  reductions  of  both  the  cavity  tones  and  the  broadband.  The  authors  attributed  the 
potential  meehanism  of  the  actuator  to  the  interaction  of  the  shed  vortices  with  the  shear  layer 
instabilities.  As  will  be  diseussed  later,  however,  there  are  additional  possible  meehanisms  noted 
by  other  investigators  that  influence  the  suppression  effectiveness  of  the  cylinder. 

At  the  same  time  in  the  mid- 1990’s,  Ahuja  and  his  colleagues  were  investigating  other 
novel  control  strategies.  For  example,  Mendoza  and  Ahuja  (1996)  studied  the  effeet  of  a  steady 
wall  jet  on  the  tone  production  meehanism,  using  a  Coanda  surfaee.  Although  no  mass  flow 
measurements  were  obtained,  upstream  boundary  layer  profiles  showed  an  increase  in  S  with 
blowing,  thereby  leading  to  the  hypothesis  that  the  suppression  was  due  to  reduced  amplification 
of  the  shear  layer  instabilities. 

Hsu  and  Ahuja  (1996)  studied  the  effeet  of  a  trailing-edge  array  of  Helmholtz  resonators 
(i.e.,  commercial  syringes)  on  cavity  noise,  and  they  obtained  some  suppression  at  lower  Mach 
numbers.  At  intermediate  Mach  numbers,  the  resonators  redueed  the  magnitude  of  the  targeted 
tone,  but  new  tones  appeared  at  other  frequencies  -  a  phenomenon  that  has  been  observed  by 
many  researchers.  At  high  Mach  numbers,  no  suppression  was  obtained,  but  the  authors 
believed  that  the  reason  was  likely  due  to  the  difficulty  in  setting  the  small  resonator  volume 
accurately.  This  study  is  noteworthy  for  its  attempt  to  eontrol  eavity  oscillations  in  the  vicinity 
of  the  aeoustie  souree  origin  near  the  trailing  edge  impingement  region.  Generally,  aetuators  are 
plaeed  at  the  leading  edge  of  the  eavity  to  leverage  the  growth  of  instabilities  in  the  shear  layer. 

Cattafesta  et  al  ( 1 997)  presented  the  use  of  a  six-element  piezoelectric  flap  array  flush 
mounted  at  the  leading  edge  of  the  eavity.  The  bandwidth  of  the  actuators  was  large  enough 
(-300  Hz)  to  provide  forcing  at  frequencies  comparable  to  that  of  the  eavity  tones.  Open-loop 
sinusoidal  foreing  at  a  sufficient  amplitude  and  appropriate  detuned  frequency  was  capable  of 
suppressing  the  eavity  tone.  Shear  layer  velocity  measurements  indicated  that  the  aetuators 
seeded  the  shear  layer  with  a  disturbance  that  was  large  enough  to  prevent  the  growth  of  the 
natural  eavity  disturbances.  However,  the  possibility  of  starving  the  growth  of  natural 
instabilities  at  high  subsonie  and  supersonic  Maeh  numbers  is  questionable,  because  of  the  large 
amount  of  mean  flow  energy  available  Ml )  for  disturbance  amplification. 

Shaw  (1998)  continued  his  study  of  leading-edge  HFTGs,  low  frequency  pulsed  fluidie 
injection,  and  oscillating  flaps.  While  various  diameter  HFTGs  mounted  at  a  fixed  height  were 
shown  to  be  effective,  the  suppression  improved  as  the  diameter  was  increased.  However,  the 
relative  height  of  the  cylinder  in  the  boundary  layer  was  not  reported,  which  is  now  known  to  be 
an  important  parameter.  Shaw  also  diseussed  two  potential  meehanisms  of  the  HFTG:  (1)  high 
frequency  foreing  due  to  shedding  and  (2)  redueed  shear  layer  growth  rates  due  to  boundary 
layer  thickening. 

Resealing  Shaw’s  pulsed  blowing  results  show  that  they  are  consistent  with  prior  results, 
since  B  values  of  a  few  pereent  were  required  to  suppress  the  tones.  However,  no  spectra  were 
reported  to  assess  the  impaet  of  blowing  on  the  broadband  noise.  Interestingly,  the  tone 
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amplitude  continued  to  decrease  as  the  pulse  frequency  of  the  injector  reached  its  upper  limit  of 
100  Hz.  Furthermore,  normal  injection  was  shown  to  be  superior  to  tangential  blowing. 

The  oscillatory  flap  frequency  in  Shaw’s  experiment  was  varied  from  0-100  Hz  and 
provided  maximum  suppression  at  5  Hz.  A  monotonic  improvement  (reduction)  in  unsteady 
pressure  level  occurred  as  the  dynamic  deflection  angle  increased  to  its  upper  limit, 
corresponding  to  a  deflection  of  order  5 .  However,  the  increase  in  S  for  a  full-scale  aircraft  led 
Shaw  to  conclude  that  this  approach  (low  frequency,  large  amplitude,  open-loop  forcing)  was  not 
feasible  for  a  full-scale  aircraft. 

In  1999,  three  new  approaches  were  reported.  Fabris  and  Williams  (1999)  used  unsteady 
bleed  (zero-net  mass-flux)  forcing  to  produce  a  broadband  actuator  capable  of  producing  a 
complex  input  disturbance  comprised  of  multiple  frequency  components.  They  demonstrated 
that  the  shear  layer  was  most  receptive  to  horizontal  or  tangential  forcing  via  shear  layer  velocity 
measurements,  in  contrast  to  the  results  of  Shaw  (1998). 

Lamp  and  Chokani  (1997)  used  a  rotary  valve  actuator  to  provide  steady  and/or 
oscillatory  blowing  upstream  of  the  cavity  leading  edge  at  a  particular  pulsing  frequency.  Their 
actuator  configuration  emphasized  three-dimensional  effects  and  showed  that  oscillatory  blowing 
can  reduce  tone  amplitudes  provided  the  forcing  frequency  is  not  a  harmonic  of  the  cavity 
resonance. 

Raman  et  al  (1999,  2004)  used  novel  miniature  fluidic  oscillators  to  suppress  cavity 
oscillations.  These  devices  were  capable  of  producing  up  to  3  kHz  oscillations  with  mass  flow 
rates  of  less  than  0.12%  of  the  main  jet  flow  and  produced  significant  tonal  reductions. 
However,  the  mass  flow  rate  and  frequency  of  oscillations  are  coupled  (albeit  in  a  predictable 
fashion).  Whether  the  steady  mass  addition  or  the  unsteady  oscillations  were  responsible  for  the 
sound  suppression  could  not  be  determined.  This  is  a  key  unresolved  issue  and,  as  emphasized 
in  Rowley  and  Williams  (2006),  independent  control  of  the  mean  and  oscillatory  components  is 
required. 

Stanck  ct  al  (2000,  2002  and  2003)  reported  on  a  series  of  larger-scalc  experiments 
conducted  in  the  United  Kingdom  over  the  past  few  years.  In  the  first  experiment  reported  in 
2000,  they  investigated  powered  resonance  tubes,  protruding  piezoceramic  driven  wedges,  a 
cylindrical  rod,  and  passive  resonance  tubes  vs.  a  conventional  spoiler.  An  interesting  result  was 
that  the  powered  resonance  tubes  demonstrated  significant  tonal  and  broadband  reduction  when 
Bc  ~  1.6%.  The  result  was  termed  a  successful  demonstration  of  high  frequency  forcing 
(defined  as  a  frequency  that  is  very  large  compared  to  that  of  the  cavity  tones). 

A  follow-on  study  reported  in  2002  investigated  powered  and  unpowered  resonance  tubes 
(in  which  the  resonator  tubes  were  blocked  to  inhibit  high  frequency  excitation),  and  microjcts 
vs.  various  other  devices.  While  the  powered  resonance  tubes  were  redesigned  to  reduce  their 
mass  flow  requirements,  optimal  suppression  still  required  B  ~  0.6%.  The  unpowered 

resonance  tubes  consistently  provided  the  best  suppression,  indicating  that  the  primary 
suppression  mechanism  of  these  devices  is  not  just  due  to  high  frequency  forcing  but  is  also 
influenced  by  the  steady  blowing  component.  The  results  also  introduced  microjet  blowing,  and 
showed  that  vertical  blowing  is  required  for  these  devices  to  be  successful  in  this  application. 

Stanek  et  al  (2003)  offered  a  new  explanation  for  the  high-frequency  forcing  effect.  The 
intrinsic  idea  was  that  high  frequency  forcing  alters  the  instability  characteristics  such  that  the 
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growth  of  large  scale  disturbances  is  inhibited  or  prevented.  They  hypothesized  that  the 
mechanism  was  a  decelerated  energy  cascade  in  contrast  to  the  findings  of  Wiltse  and  Glezer 
(1998). 

In  2003,  Stanek  et  al.  reported  various  aspects  of  the  cylindrical  rod  in  crossflow.  They 
studied  the  vertical  position  of  the  rod  H j8  in  the  boundary  layer,  its  relative  size  d/S , 
installation  issues,  and  end  conditions.  They  recommended  an  optimal  location  as  centered  at  the 
edge  of  the  boundary  layer  and  an  optimal  size  of  d/S  =  2/3.  They  argued  that  their  results 
conclusively  demonstrate  that  the  suppression  is  due  to  high  frequency  forcing  via  vortex 
shedding  from  the  cylinder.  Additional  arguments  have  also  been  provided  that  diseuss  the  loeal 
stabilization  of  the  shear  layer  via  high  frequency  forcing.  While  the  cylinder  elearly  affects  the 
mean  flow  and  its  stability  characteristics,  there  are  other  important  faetors  that  cannot  be 
ignored,  including  experimental  evidence  presented  by  Ukeiley  et  al  (2004)  and  the  numerical 
simulations  of  Arunajatesan  et  al  (2003)  which  indicate  that  the  cylinder  ean  also  lift  the  shear 
layer  and  cause  the  impingement  region  to  be  altered.  If  the  shear  layer  impingement  location  is 
indeed  altered,  then  the  source  strength  is  presumably  affected. 

It  is  clear  that  these  studies  provide  insufficient  information  to  sort  out  these  different 
physical  mechanisms.  To  do  so  requires  confirmation  in  the  form  of  detailed  experiments, 
analysis,  and  validated  simulations  to  determine  the  mean  velocity  profile  and  subsequent  shear 
layer  instability  characteristics  for  various  high  frequency  devices. 

There  are  a  few  other  studies  involving  steady  and/or  pulsed  blowing  that  have  provided 
physical  insight  or  have  shown  promising  results.  Bueno  et  al  (2002)  used  an  array  of  six  fast- 
response  (~3  ms)  miniature  jets  mounted  upstream  of  the  leading  edge  to  study  the  effeets  of 
normal  injection  on  a  Maeh  2  cavity  flow.  They  used  instantaneous  and  ensemble-averaged 
pressure  time  histories  and  cross  correlations  to  study  the  effects  of  single  short  and  long  cyclical 
pulses  (50%  duty  eyele),  the  latter  with  relatively  low  foreing  frequencies  (50  or  80  Hz) 
compared  to  that  of  the  tones.  They  compared  their  pulsed  results  to  steady  blowing  with  B  = 
0.28%,  0.24%,  0.18%,  and  0.16%  at  L/D=  5,  6,  8,  and  9,  respectively,  and  concluded  that 
continuous  mass  injection  is  more  effective  for  suppression  than  low  frequency  pulsed  blowing. 

Ukeiley  et  al  (2003)  used  an  array  of  eight  powered  whistles  mounted  in  the  forward 
cavity  wall  as  flow  eontrol  aetuators.  These  deviees  essentially  produee  a  high  frequency  tonal 
oseillation  superimposed  on  a  steady  jet.  The  jet  is  directed  in  the  downstream  direction  but  has 
a  slight  vertical  velocity  component.  The  authors  studied  the  novel  use  of  different  injection 
gases  (heated  air,  nitrogen,  and  helium)  with  and  without  the  high  frequency  “whistle” 
component.  The  best  sound  suppression  results  were  obtained  using  steady  helium  blowing  (no 
high  frequency  component)  with  very  low  B  =  0.09%.  The  suppression  mechanism  requires 

further  study,  but  sample  Particle  Image  Veloeimetry  (PIV)  images  and  cross  correlations  of 
pressure-time  histories  suggest  that  the  injection  alters  the  impingement  region  and  disrupts  the 
acoustic  feedback  loop.  Their  results  also  highlight  the  need  to  rigorously  study  high  frequency 
foreing  effeets  isolated  from  the  mass  flow  injection  from  the  actuator. 

Zhuang  et  al  (2003)  investigated  the  use  of  a  vertically  directed  mierojet  array  mounted 
upstream  of  the  cavity  leading  edge.  The  mierojets  had  a  400  “m  diameter  and  produced  sonie 
jets  that  interact  with  the  upstream  boundary  layer.  The  authors  show  how,  at  Maeh  2,  the 
mierojets  alter  the  cavity  shear  layer  thickness  and  the  reeeptivity  process,  introduce  streamwise 
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vortieity,  and  alter  the  shear  layer  trajectory  and  the  resulting  impingement  region.  Significant 
tonal  and  broadband  suppression  levels  were  achieved  with  Bc  as  low  as  0.15%,  which  is 
significantly  lower  than  the  mass  injection  required  in  other  studies  using  steady  blowing. 
Higher  levels  of  Bc  produced  no  significant  improvement. 

Collectively,  the  blowing  results  described  above  indicate  that  manually  optimized  steady 
blowing  configurations  with  B  <  0.2%  can  be  effective  suppression  devices.  At  subsonie 

speeds  the  primary  mechanisms  appear  to  be  an  alteration  of  the  shear  layer  stability 
characteristics,  the  introduction  of  streamwise  vortieity,  and  shear  layer  impingement  location. 
While  these  arc  also  important  at  supersonic  speeds,  shock  wave/boundary  layer  interactions  at 
the  upstream  cavity  edge  and  the  ensuing  shear  layer  trajectory  alteration  appear  to  be  additional 
factors  that  should  be  considered. 

It  is  interesting  that  when  all  of  the  available  blowing  data  is  expressed  using  the  blowing 
coefficient  definition  of  Vakili  and  Gauthier  (1994)  one  finds  that  the  evolution  of  steady 
blowing  techniques  has  reduced  the  effective  (not  necessarily  optimal)  Bc  from  0(10%)  by  two 
orders  of  magnitude  down  to  0(0. 1%). 

High  frequency  excitation,  whether  it  is  passive  or  active,  appears  promising  for  both 
tonal  and  broadband  suppression.  However,  the  responsible  mechanisms  require  further  study. 
There  is  ample  evidence  that  high  frequency  foreing  alters  the  mean  flow.  As  a  result,  the  shear 
layer  stability  characteristics  are  altered  and,  in  some  eases,  the  trajectory  of  the  shear  layer  is 
modified.  When  the  impingement  location  is  altered,  the  strength  of  the  aeoustie  souree  is 
redueed  and  the  broadband  noise  level  decreases.  To  date  closed-loop  control  produces 
comparatively  little  ehange  in  the  mean  flow  properties  and,  as  such,  has  only  been  shown  to  be 
effective  for  tonal  control. 


1.2.3  Closed  Loop  Adaptive  Control  Approaches 

Aetivc  feedback  flow  eontrol  has  been  applied  to  the  flow-induced  cavity  oscillation 
problems  over  approximately  the  past  20  years.  The  closed-loop  eontrol  approaches  have 
advantages  of  redueed  energy  consumption  (Cattafesta  et  al.  1997),  no  additional  drag  penalty, 
and  robustness  to  parameter  changes  and  modeling  uncertainties.  In  general,  elosed-loop  flow 
eontrol  measures  and  feeds  back  pressure  fluctuations  at  the  surface  of  the  cavity  walls  (or  floor) 
to  an  actuator  at  the  cavity  leading  edge  to  suppress  the  cavity  oscillations  in  a  elosed-loop 
fashion. 

In  general,  past  active  control  strategies  have  taken  one  of  two  approaches  for  the 
purpose  of  reducing  cavity  resonance.  First,  they  can  thicken  the  boundary  layer  in  order  to 
reduce  the  growth  of  the  instabilities  in  the  shear  layer.  Alternatively,  they  ean  be  used  to  break 
the  internal  feedback  loop  of  the  cavity  dynamics.  Most  closed-loop  schemes  exploit  the  latter 
approach.  Early  closed-loop  control  applications  used  manual  tuning  of  the  gain  and  delay  of 
simple  feedback  loops  to  suppress  resonance  (Gharib  et  al.  1987;  Williams  et  al.  2000a, b). 
Mongeau  et  al.  (1998)  and  Kook  et  al.  (2002)  used  an  active  spoiler  driven  at  the  leading  edge 
and  a  loop-shaping  algorithm  to  obtain  significant  attenuation  with  small  actuation  effort. 
Debiasi  et  al.  (2003,  2004)  and  Samimy  et  al.  (2003)  proposed  a  simple  logic-based  controller  for 


closed-loop  cavity  flow  control.  Low-order  model-based  controllers  with  different  bandwidths, 
gains  and  time  delays  have  also  been  designed  and  implemented  (Rowley  et  al.  2002,  2003, 
Williams  et  al.  2002,  Micheau  et  al.  2004,  Debiasi  et  al.  2004).  Linear  optimal  controllers 
(Cattafesta  et  al.  1997,  Cabell  et  al.  2002,  Debiasi  et  al.  2004,  Samimy  et  al.  2004,  Caraballo  et 
al.  2005)  have  been  successfully  designed  for  operation  at  a  single  flow  condition.  These  models 
are  all  based  on  reduced-order  system  models,  and  most  of  these  controller  design  methods  are 
based  on  model  forms  of  the  frequency  response  function,  rational  diserete/eontinuous  transfer 
function,  or  state-space  form.  However,  the  coefficients  of  these  model  forms  are  assumed  to  be 
constant,  and  this  assumption  requires  that  the  system  is  time  invariant  or  at  least  a  quasi-static 
system  with  a  fixed  Mach  number. 

Although  the  physical  models  of  flow-induced  cavity  oscillations  have  been  explored 
extensively,  they  are  not  convenient  for  control  realization.  This  is  because  these  models  are 
highly  dependent  on  the  accuracy  of  the  estimated  internal  states  of  the  cavity  system.  In 
addition,  cavity  flow  is  known  to  be  quite  sensitive  to  slight  changes  in  flow  parameters.  So  a 
small  change  in  Mach  number  can  deteriorate  the  performance  of  a  single-point  designed 
controller  (Rowley  and  Williams  2003).  Therefore,  adaptive  control  is  certainly  a  reasonable 
approach  to  consider  for  reducing  oscillations  in  the  flow  past  a  cavity.  Adaptive  control 
methodology  combines  a  general  control  strategy  and  system  identification  (ID)  algorithms. 
This  method  is  thus  potentially  able  to  adapt  to  the  changes  of  the  cavity  dimension  and  flow 
conditions.  It  updates  the  controller  parameters  for  optimum  performance  automatically. 

Two  distinct  loops  are  typically  observed  in  the  controller  using  this  method.  The  outer 
loop  is  a  standard  feedback  control  system  comprised  of  the  process  block  and  the  controller 
block.  The  controller  operates  at  a  sample  rate  that  is  suitable  for  the  discrete  process  under 
control.  The  inner  loop  consists  of  a  parameter  estimator  block  and  a  controller  design  block. 
An  ID  algorithm  and  a  specified  cost  function  are  then  used  to  design  a  controller  that  will 
minimize  the  output.  The  steps  for  real-time  flow  control  include:  (i)  Use  a  broadband  system 
ID  input  from  the  actuator(s)  and  the  measured  pressure  fluctuation  output(s)  on  the  walls  of  the 
cavity  to  estimate  the  system  (plant  and  disturbance)  parameters,  (ii)  Design  a  controller  based 
on  the  estimated  system  parameters,  (iii)  Control  the  whole  system  to  minimize  the  effects  of  the 
disturbance,  measured  noise,  and  the  uncertainties  in  the  plant. 

Based  on  this  adaptive  control  methodology,  some  adaptive  algorithms  adjust  the 
controller  design  parameters  to  track  dynamic  changes  in  the  system.  However,  only  a  few 
researchers  have  demonstrated  the  on-line  adaptive  closed-loop  control  of  flow-induced  cavity 
oscillations.  Cattafesta  et  al.  (1999  a,  b)  applied  an  adaptive  disturbance  rejection  algorithm, 
which  was  based  upon  the  ARMARKOV/Toeplitz  models  (Akers  and  Bernstein  1997; 
Venugopal  and  Bernstein  2000,  2001),  to  identify  and  control  a  cavity  flow  at  Mach  0.74  and 
achieved  10  dB  suppression  of  a  single  Rossiter  model.  Other  modes  in  the  cavity  spectrum 
were  unaffected.  Insufficient  actuator  bandwidth  and  authority  limited  the  control  performance 
to  a  single  mode.  Williams  and  Morrow  (2001)  applied  an  adaptive  filtcrcd-X  LMS  algorithm  to 
the  cavity  problem  and  demonstrated  multiple  cavity  tone  suppression  at  Mach  number  up  to 
0.48.  However,  this  was  accompanied  by  simultaneous  amplification  of  other  cavity  tones. 
Numerical  simulations  using  the  least  mean  squares  (LMS)  algorithm  were  shown  by  Kestens 
and  Nieoud  (1998)  to  minimize  the  output  of  a  single  error  sensor.  The  reduction  was  associated 
with  a  single  Rossiter  mode,  but  only  within  a  small  spatial  region  around  the  error  sensor. 
Kegcrise  et  al.  (2002)  implemented  adaptive  system  ID  algorithms  in  an  experimental  cavity 
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flow  at  a  single  Mach  number  of  0.275.  They  also  summarized  the  typical  finite-impulse 
response  (FIR)  and  infinite-impulse  response  (IIR)  based  system  ID  algorithms.  They  concluded 
that  the  FIR  filters  used  to  represent  the  flow-induced  cavity  process  were  unsuitable.  On  the 
other  hand,  IIR  models  were  able  to  model  the  dynamics  of  the  cavity  system.  LMS  adaptive 
algorithm  was  more  suitable  for  real-time  control  than  the  recursive-least  square  (RLS)  adaptive 
algorithm  due  to  its  reduced  computational  complexity.  Recently,  more  advanced  controllers, 
such  as  direct  and  indirect  synthesis  of  the  neural  architectures  for  both  system  ID  and  control 
(Efe  et  al.  2005)  and  the  generalized  predictive  control  (GPC)  algorithm  (Kegerisc  et  al.  2004), 
have  been  implemented  on  the  cavity  problems. 

From  a  physical  point  of  view,  the  closed-loop  controllers  have  no  effect  on  the  mean 
velocity  profile  (Cattafesta  et  al.  1997).  However,  they  significantly  affect  streamwisc  velocity 
fluctuation  profiles.  This  control  effect  eliminates  the  strength  of  the  pressure  fluctuations 
related  to  flow  impingement  on  the  trailing  edge  of  the  cavity.  Although  closed-loop  control  has 
provided  promising  results,  the  peaking  (i.e.,  generation  of  new  oscillation  frequencies),  peak 
splitting  (i.e.,  a  controlled  peak  splits  into  two  sidebands)  and  mode  switching  phenomena  (i.e., 
non-linear  interaction  between  two  different  Rossiter  frequencies)  often  appear  in  active  closed- 
loop  control  experiments  (Cattafesta  et  al.  1997,  1999  b;  Williams  et  al.  2000;  Rowley  et  al.  2002 
b,  2003;  Cabell  et  al.  2002;  Kegerisc  et  al.  2002,  2004a). 

Explanations  of  these  phenomena  are  provided  by  Rowley  et  al.  (2002b,  2006),  Banaszuk 
et  al.  (1999),  Hong  and  Bernstein  (1998),  and  Kegerisc  et  al.  (2004).  Rowley  et  al.  (2002b, 
2003)  concluded  that  if  the  viewpoint  of  a  linear  model  was  correct,  a  closed-loop  controller 
could  not  reduce  the  amplitude  of  oscillations  at  all  frequencies  as  a  consequence  of  the  Bode 
integral  constraint.  Banaszuk  et  al.  (1999)  gave  explanations  of  the  peak-splitting  phenomenon. 
They  claimed  that  the  peak  splitting  effect  was  caused  by  a  large  delay  and  a  relatively  low 
damping  coefficient  of  the  open-loop  plant.  Cabell  et  al.  (2002)  explained  these  phenomena  by 
the  combination  of  inaccuracies  in  the  identified  plant  model,  high  gain  controllers,  large  time 
delays  and  uncertainty  in  system  dynamics.  In  addition,  narrow-bandwidth  actuators  and 
controllers  may  also  lead  to  a  peak-splitting  phenomenon  (Rowley  et  al.  2006). 

Hong  and  Bernstein  defined  the  closed-loop  system  disturbance  amplification  (peaking) 
phenomenon  as  spillover.  They  illustrated  that  the  spillover  problem  was  caused  by  the 
collocation  of  disturbance  source  and  control  signal  or  the  collocation  of  the  performance  and 
measurement  sensors.  For  this  reason,  the  reduction  of  broadband  pressure  oscillations  was  not 
possible  if  the  control  input  was  collocated  with  the  disturbance  signal  at  the  leading  edge  of 
cavity.  Therefore,  Kegerise  et  al.  (2004)  suggested  a  zero  spillover  controller  which  utilized 
actuators  at  both  the  leading  and  trailing  edges  of  the  cavity  for  closed-loop  flow  control. 


Objectives 

There  were  several  objectives  of  this  study  which  were; 

•  Develop  methodology  that  will  allow  for  the  study  of  time  dependent  features. 
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•  Formulate  an  approach  that  can  utilize  time  dependent  velocity  and  density  field 
information  to  investigate  the  sources  of  the  large  fluctuating  pressures  observed  in 
eavities. 

•  Develop  a  tool  that  ean  be  used  for  quantitatively 
be  used  in  aeeordanee  with  the  estimation  procedure. 

•  Refine  and  apply  a  closed  loop  flow  eontrol 
effectiveness 

•  Conduct  studies  to  examine  the  differences  in  the  flow  structure  between  natural 
and  effectively  controlled  flows 


asses  the  density  field  that  ean 
strategy  to  demonstrates  the 


1.4  Layout  of  this  Document 

In  this  document  we  will  discuss  our  aeeomplishments  during  the  course  of  this  AFOSR 
grant  all  of  whieh  are  pieees  that  will  allow  for  us  to  efficiently  reduee  the  surfaee  pressure 
fluctuations  in  open  eavities  in  supersonie  flows.  The  remaining  parts  of  the  document  contain  a 
ehapter  on  each  one  of  these  pieees.  Chapter  2  presents  the  methodology  and  application  for  the 
estimation  of  the  time  dependent  vcloeity  field  based  on  surfaee  pressure.  Additionally  we  will 
present  a  method  for  method  of  integrating  the  velocity  field  with  Poisson’s  equation  to  calculate 
the  different  components  of  the  sources  of  the  fluctuating  surface  pressure.  In  Chapter  3  data 
acquired  in  both  a  controlled  and  uncontrolled  Mach  1.5  flow  over  an  open  cavity  is  analyzed 
w  ith  the  Proper  Orthogonal  Decomposition  to  evaluate  similarities  in  the  flow'  structure.  Chapter 
4  demonstrates  the  “state  of  the  art”  in  adaptive  flow  eontrol  for  surfaee  pressure  reductions  in 
flow  over  a  eavity  along  with  a  zero  net  mass  flux  aetuator.  In  Chapter  5  we  present  the 
development  and  application  of  optical  defleetometry  whieh  in  the  future  ean  be  used  in 
conjunction  with  synchronous  surface  pressure  fluctuation  measurements  for  the  time  dependent 
flow  field  estimation  procedure. 
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2  Estimation  of  Time  Dependent  Phenomena 


2.1  Methodology 

Redueed  order  descriptions  are  needed  to  describe  eomplex  flows,  like  that  over  a  cavity, 
so  that  the  underlying  physical  mcehanisms  can  be  interpreted.  Over  the  past  deeade, 
methodologies  for  estimating  low  dimensional  descriptions  of  a  time-resolved  flow  field  using 
boundary  information  have  proven  to  be  instrumental  in  fluid  dynamics.  Some  of  the  more 
promising  methodologies  arc  rooted  in  large  scale  turbulent  structure  detection  techniques  sueh 
as  Stochastie  Estimation  (Adrian  1975)  and  the  Proper  Orthogonal  Decomposition  (Lumlcy 
1967). 

The  Proper  Orthogonal  Decomposition  (POD)  is  a  mathematically  unbiased  technique  for 
decomposing  a  flow  field  into  an  orthonormal  basis  set  whieh  is  optimal  in  a  mean  square  sense. 
This  basis  set  is  found  from  the  solution  of  the  following  integral  eigenvalue  problem, 

j 3  (x,x)(p"  (x')cix'  =  TV,"  (*)  (2. 1 ) 

where  T  is  the  eignevalue  and  represents  the  weight  of  each  given  mode  and  <J)  are  the  orthogonal 
eigenfunetions  whieh  represent  the  spatial  distribution  of  the  basis  set.  The  kernel  R,,  is  the  two- 
point  spatial-correlation  tensor  and  can  be  formulated  for  sealar  or  veetor  flow  properties  sueh  as 
the  density  and  velocity  fields  that  will  be  discussed  below.  Hilbert-Schmidt  (see  Holmes  ct  al. 
1996)  theory  specifies  that  there  is  a  denumerable  set  of  the  eigenmodes  that  are  orthonormal, 
and  that  the  instantaneous  field,  in  this  example  velocity,  can  be  represented  as  a  linear 
combination  of  the  product  of  time  dependent  coefficients  with  spatial  basis  functions 

=  (2.2) 

n 

where  the  expansion  coefficients  are  ealeulated  from 

a"  =  {xj)<j>"  (x)dx .  (2.3) 

n 

In  order  to  obtain  the  time  dependence  of  the  expansion  coefficients,  the  approach 
requires  that  the  time-resolved  instantaneous  field  is  known  at  all  spatial  locations 
simultaneously  to  solve  the  integral  in  equation  2.3.  This  condition  is  not  usually  met  in 
experiments,  and  one  must  use  alternative  methods,  such  as  the  Stochastic  Estimation  or  low 
dimensional  models  based  on  Galerkin  projections  (Ukciley  et  al.  2001)  to  determine  the  time 
dependence  of  the  expansion  coefficients. 

The  use  of  Stochastic  Estimation  combined  with  the  POD  has  been  termed 
“ modified  Stochastic  Estimation”  (mQSE)  in  the  literature  and  has  been  used  extensively  for 
studying  the  time  dependent  features  of  subsonic  cavity  flow  (Ukeiley  and  Murray  2005  and 
Murray  and  Ukeiley  2007).  Perhaps  one  of  the  most  important  and  novel  features  of  the 
application  in  those  studies  was  the  use  of  the  surfaee  pressure  as  the  estimation  condition.  In 
particular,  sinee  time-resolved  pressure  measurements  are  readily  obtained,  this  results  in  a  time- 
resolved  estimate  of  the  low-dimensional  velocity  field.  The  estimation  can  be  written  as 
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(2.4) 


a'(/)  =  (a"(/)lA(')) 

which  states  that  the  POD  expansion  coefficient  is  estimated  given  statistical  information  about 
its  relationship  to  the  surface  pressure  at  several  locations  k.  Just  as  with  the  formulation  of  the 
conventional  stochastic  estimation  procedure,  the  estimated  POD  expansion  coefficient  can  be 
written  in  a  series  expansion  as 

5"(t)  =  A,pt(t)  +  B„lmPl(t)pJt)  (2.5) 


which  has  been  truncated  at  the  quadratic  term  (see  Murray  and  Ukeiley  2003  for  a  validation). 
Minimizing  the  mean-square  error  between  the  estimated  and  true  expansion  coefficient  results 
in  a  matrix  equation 


M  «[/»]>]• 


(2.6) 


where  the  matrices  represent  the  various  correlations  between  the  surface  pressure  with  itself  and 
the  surface  pressure  with  the  POD  expansion  coefficients  and  their  explicit  definitions  can  be 
found  in  Murray  and  Ukeiley  (2007).  Finally,  the  time  dependent  estimated  velocity  field  can  be 
reconstructed  solving  equation  2  with  the  solutions  from  equation  5. 


Figure  1  summarizes  the  estimation  procedure  that  is  used  to  obtain  the  time 
dependent  low-dimensional  description  of  the  cavity  flow  field.  The  top  of  this  figure 
demonstrates  that  one  starts  from  time  resolved  surface  pressure  data  and  independent  snapshots 
of  either  the  velocity  or  density  fields.  Then,  through  the  equations  described  above,  one 
estimates  a  low-dimensional  representation  of  these  fields  as  depicted  by  the  dashed  red  lines  in 
the  bottom  part  of  the  figure. 
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Figure  1:  Flow  Chart  of  Time  Dependent  Estimation  Procedure. 


It  can  be  noted  that  in  the  formulation  above  we  are  using  a  static  estimator  of  the  flow 
field.  However  in  the  applications  discussed  below,  the  array  of  surface  pressure  sensors  spans 
the  whole  domain  of  interest,  which  overcomes  the  need  to  account  for  convection  in  the  flow 
field.  However,  we  will  present  work  below  that  uses  a  dynamic  estimator  of  the  surface 
pressure  to  incorporate  memory  effects. 
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2.2  Application 

There  were  two  sets  of  pressure  data  utilized  in  performing  the  mQSE.  The  first  set, 
Pk(te),  arc  the  pressure  values  associated  with  the  0=1  to  M  snapshots  of  the  flow  field.  These 
values  arc  used  to  calculate  the  correlations  necessary  to  compute  the  estimation  coefficients.  If 
Pk(to)  is  then  reused  to  estimate  the  velocity  field  at  every  0,  the  estimation  effectively  attempts 
to  reconstruct  the  P1V  data.  This  allows  the  estimated  flow  field  and  single-point  statistics  to  be 
compared  with  the  original  data.  The  second  set  of  pressure  data,  Pk(te),  represent  0=1  to  N 
time-resolved  pressure  measurements  sampled  at  90  kHz.  This  data  was  used  to  estimate  the  time 
dependence  of  the  flow  using  the  mQSE.  The  results  of  this  estimation  arc  presented  for  several 
free  stream  Mach  numbers.  For  all  the  estimates,  the  mQSE  is  applied  only  to  the  fluctuating 
quantities. 

Evaluation  of  the  Velocity  Estimation 

Based  on  the  results  from  the  Gappy  POD  to  the  experimental  data  set  used  here  (see 
Murray  and  Ukeiley,  2007),  only  the  first  8  POD  modes  were  used  in  the  mQSE  to  estimate  the 
velocity  field.  The  fraction  of  the  total  kinetic  energy  represented  by  these  first  8  POD  modes 
can  be  calculated  by  summing  over  the  respective  eigenvalues.  At  each  Mach  number,  the  first  8 
POD  modes  account  for  approximately  50%  of  the  total  energy.  Therefore,  in  order  to  match  the 
kinetic  energy  present  in  the  original  data,  the  estimated  velocity  fluctuations  were  scaled  up  by 
the  reciprocal  of  the  include  energy  fraction.  For  example,  if  8  modes  represented  49%  of  the 
total  energy,  the  estimated  velocity  fluctuations  were  scaled  up  by  a  factor  of  2.04. 


m 


Figure  2:  Comparison  of  the  turbulence  intensities  between  the  PIV  data  and  the  mQSE 
estimation  of  the  flow  field  for  Maeh  0.58:  (a),  (e),  &  (e)  are  streamwise  turbulence  intensity  and 
(b),  (d),  &  (f)  are  wall-normal  turbulence  intensity,  (a)  &  (b)  are  from  the  PIV  data,  (c)  &  (d) 
were  calculated  from  an  8  mode  POD  reconstruction  of  the  PIV  data,  and  (e)  &  (f)  are  from  the 
mQSE  estimation. 
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Figure  2  compares  the  streamwise  and  wall-normal  turbulence  intensities  resulting  from  the 
mQSE  estimation  with  the  statistics  of  the  original  PIV  data.  Plots  (a)  and  (b)  show  the 
turbulence  intensities  for  the  PIV  data.  In  plots  (c)  and  (d),  the  statistics  wrere  calculated  for  the 
PIV  data  after  first  reconstructing  the  snapshots  using  the  first  8  POD  modes.  Plots  (e)  and  (f) 
show  the  statistics  calculated  from  the  mQSE  estimation  of  the  PIV  data.  The  contour  levels 
show  general  agreement  in  both  structure  and  magnitude  between  the  original  data  and  the 
mQSE  estimation.  Similar  results  were  obtained  for  the  other  Mach  numbers.  A  qualitative 
evaluation  of  the  estimate  is  presented  in  figure  8  where  a  PIV  snapshot  taken  from  the  Mach 
0.58  data  set  is  estimated  using  the  associated  pressure  measurements  plotted  at  the  bottom  of  the 
figure.  Plot  (a)  shows  the  original  PIV  data,  and  plot  (b)  shows  an  8  mode  POD  reconstruction  of 
the  same  snapshot.  In  plot  (c)  this  snapshot  has  been  estimated  with  the  mQSE  using  the  pressure 
distribution  in  plot  (d).  In  all  the  plots  of  the  flow  field,  the  contours  which  outline  the  spanwise 
vorticity  have  the  same  magnitudes  and  the  vectors  arc  plotted  using  the  same  reference  length. 
The  validity  of  the  chosen  scaling  parameter  (based  on  the  percentage  of  energy  included  in  the 
first  8  POD  modes)  is  demonstrated  because  the  strength  of  the  spanwise  vorticity  in  (c) 
compares  more  favorably  with  the  original  PIV  data  in  (a)  than  the  8  mode  POD  reconstruction 
shown  in  (b).  However,  the  structure  of  the  vorticity  field  in  (c)  is  very  similar  to  that  in  (b) 
which  demonstrates  the  filtering  that  results  from  using  only  the  first  8  POD  modes  in  the 
mQSE. 

Time  Dependent  Estimation  of  the  Flow  Field 

The  results  presented  in  this  section  demonstrate  the  temporal  dynamics  of  the  flow  field 
estimated  using  the  mQSE.  The  main  focus  here  is  the  comparison  between  the  flow  field 
dynamics  and  the  instantaneous  surface  pressure  loads  which  are  shown  along  with  each 
realization  of  the  flow. 

Each  figure  in  this  section  shows  nine  consecutive  realizations  of  the  flow  field.  They  arc 
arranged  such  that  time  proceeds  from  top-to-bottom  and  left-to-right  in  each  figure  starting  at  t« 
and  following  at  the  noted  intervals.  The  wall-pressure  distribution  is  plotted  directly  below  each 
image  of  the  flow  field  using  a  bar-chart  with  a  connecting  spline  curve.  The  discussion  that 
follows  also  takes  into  account  observations  made  by  examining  extended  time-series 
animations. 

Figure  3  through  Figure  7  show  the  temporal  evolution  of  the  velocity  field  estimated 
using  the  mQSE.  Contours  of  clockwise  rotating  spanwise  vorticity  arc  also  shown  in  the  plots. 
Figure  3  shows  the  mQSE  results  for  the  Mach  0.19  flow.  Examination  of  the  estimated  time- 
series  does  not  demonstrate  a  clear  downstream  propagation  of  structures  in  the  shear  layer  (a 
downstream  motion  could  only  be  noticed  sporadically).  Additionally,  the  surface  pressure  loads 
do  not  exhibit  a  clear  wave-like  structure.  However,  it  is  evident  that  increases  in  the  wall 
pressure  on  the  cavity  floor  occur  below  regions  where  the  shear-layer  streamlines  are  curved 
away  from  the  cavity  while  decreases  in  the  pressure  accompany  curvatures  toward  the  cavity. 
The  stream-traces,  which  were  plotted  using  Tccplot’s  built-in  utility,  demonstrate  the  prevailing 
motion  that  occurs  in  the  flow:  The  circulation  region  inside  the  cavity  oscillates  back  and  forth 
without  any  clear  periodicity.  Likewise,  the  aft-wall  pressure  loads,  which  are  larger  than  at  any 
other  measurement  location,  were  found  to  oscillate  from  large  negative  to  large  positive  without 
a  well-defined  regularity.  However,  it  did  appear  that  the  aft  wall  pressure  was  related  to  the 
shear  layer  streamlines  either  stagnating  on  or  lifting  over  the  aft  wall.  Figure  4  shows  very 
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similar  results  for  the  Mach  0.29  flow.  At  this  Maeh  number  the  shear  layer  still  shows 
intermittent  patterns  of  downstream  propagating  structures.  In  contrast  to  the  Mach  0.19  flow, 
this  intermitteney  is  increased  which  explains  the  presence,  while  small  in  magnitude,  of  Rossiter 
tones  in  the  measured  spectra.  This  demonstrates  the  onset  of  well  organized  cavity  flow 
oscillations  that  will  become  well  defined  at  Mach  0.39.  These  results  suggest  that  fully 
developed  oscillations  become  sustained  for  the  current  configuration  around  Maeh  0.3. 


Figure  3:  mQSE  estimation  of  the  Mach  0.19  cavity  flow.  At  =  44ps.  Contours  outline 
areas  of  large  negative  spanwise  vorticity.  Streamtraccs  were  plotted  using  Tecplot’s  streamtrace 
utility.  Pressure  values  are  in  psi  and  plotted  at  their  respective  x/D  location. 
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Figure  4 :  Estimation  of  Mach  0.29  cavity  flow:  At  =  33(is. 


Figure  5:  Estimation  of  Mach  0.39  cavity  flow:  At  =  22|is. 
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Figure  6:  Estimation  of  Mach  0.58  cavity  flow:  At  =  1 1  fj.s. 


Figure  7:  Estimation  of  Mach  0.73  cavity  flow:  At=  1 1  (as. 
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Figure  5,  Figure  6  and  Figure  7  show  the  results  for  Maeh  0.39,  0.58,  and  0.73 
respectively.  Here,  the  stream-traees  have  been  removed  beeause  the  formation  and  downstream 
propagation  of  spanwise  vortex  structures  is  elearly  evident  as  shown  by  the  eontour  lines.  In 
eaeh  figure,  a  structure  is  labeled  as  A  and  followed  as  it  moves  downstream  through  the  shear 
layer.  The  location  of  the  structure  is  always  aeeompanied  by  a  negative  pressure  fluctuation  on 
the  eavity  floor  below  it.  The  dominance  of  a  particular  Rossiter  mode  is  observed  in  these 
figures  by  noticing  the  number  of  structures  that  oeeupy  the  shear  layer  at  any  instant.  At  Maeh 
0.39,  there  are  generally  two  structures  with  a  third  either  forming  or  passing  over  the  aft  wall. 
This  is  consistent  with  Rossiter  mode  3  whieh  was  shown  to  be  dominant  at  this  Maeh  number. 
Similarly,  at  Maeh  0.58  and  0.73,  there  is  generally  only  one  structure  with  another  being  either 
formed  or  interacting  with  the  aft  wall.  This  is  consistent  with  Rossiter  mode  2  whieh  is  the 
dominant  mode  for  these  Maeh  numbers.  The  estimated  temporal  evolution  of  the  higher  Maeh 
numbers  demonstrates  the  relationship  between  the  vortex  structures  and  the  aft-wall  pressure 
loads.  Positive  pressure  fluctuations  oeeur  just  after  a  vortex  structure  has  passed  over  the  aft 
wall  and  the  shear  layer  is  then  pulled  into  the  eavity  by  the  circulation  of  the  approaehing 
vortex.  This  eauses  the  shear  layer  to  impaet  the  aft  wall  and  stagnate  there.  Then,  as  the 
approaehing  vortex  reaehes  the  aft  wall  it  brings  the  low  pressure  at  its  eore  along  with  it  while 
also  deflecting  the  shear  layer  over  the  aft-wall  resulting  in  a  decrease  in  the  pressure. 


2.3  Source  Identification 


2.3.1  Methodology 


An  explicit  equation  for  the  pressure  in  terms  of  the  velocity  ean  be  written  for  the 
incompressible  flow  by  taking  the  divergence  of  the  momentum  equation  and  eliminating  terms 
using  the  continuity  requirement.  The  result  is  Poisson’s  equation: 


V2P  =  -P 


d2 


dxi  dxj 


-U,u,. 


(2.7) 


Taking  eues  from  the  work  by  Kraiehnan  (1956)  the  above  equation  ean  be  integrated  to  give  the 
instantaneous  pressure,  P(x,  t),  at  a  point,  x,  on  a  solid  surfaee:, 


p(rt\  =  2^-f  y 

[,)  4*  yv2>0  dxidxj  |f-r 


(2.8) 


where  the  integration  is  earned  out  over  the  positive  half-plane  (y^  >  0)  where  the  flow  exists  and 
a  multiplicative  faetor  of  2  is  included  so  that  the  flow  (souree  region)  is  effectively  mirrored 
aeross  the  bounding  surfaee. 

Setting  x  =  0,  defining  r  =  |x-y|,  and  performing  a  double  partial  integration  on  the 
double  divergence  of  the  velocity  produet  yields 
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(2.9) 


~  r2Sjj 


“V— 

Qa 

The  term  Q,j  is  the  classic  form  for  a  traceless  tensor  and  drops  off  rapidly  with  l/r\  If  the  cross- 
stream  velocity,  U3,  is  zero,  five  of  the  nine  terms  in  Q,jR,j  arc  eliminated  which  simplifies  the 
integral: 


)Ui(y%t)Uj{y,t). 
S - - - - 


With  U3  =  0,  the  quadrupole  becomes  aligned  with  the  streamwise  direction  and  inclined  by  45 
degrees  to  the  horizontal.  However,  this  does  not  completely  remove  the  dependence  of  the 
calculated  pressure  on  the  cross-stream  coordinate,  y3,  because  it  is  still  included  in  the  distance 
r2.  Figure  8  shows  Q,j  plotted  relative  to  the  physical  dimensions  of  the  cavity  showing  the 
orientation  and  extent  of  the  quadrupole  for  a  calculation  of  the  surface  pressure  at  the  center  of 
the  cavity  floor.  The  figure  demonstrates  how  fast  the  magnitude  drops  of  w  ith  1/r5  and  shows 
that  the  effects  of  the  end  walls  w  ill  be  negligible  for  locations  near  the  center  of  the  cavity. 


9 
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Figure  8:  Depiction  of  the  quadrupole  Qij  plotted  relative  to  the  physical  dimensions  of 
the  rectangular  cavity.  Flow  would  be  in  the  positive  streamwise  direction.  The  position  of  the 
quadrupole  in  the  cavity  is  appropriate  for  calculating  the  surface  pressure  at  the  center  of  the 
cavity:  (a)  shows  the  full  3-dimensional  extentand  orientation  of  the  quadrupole  in  the  cavity  for 
U3  =  0,  (b)  shows  a  2-D  slice  along  the  cavity  centerline. 
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Decomposing  the  pressure  and  velocity  into  mean  and  fluctuating  components  allows 
equation  (2.7)  to  be  written  for  the  fluctuating  pressure, 


V  p—  pQx.dx  (f  1  U->  +  ^  jU,,+  UtUj  lit  uj) '  (2.11) 

1  2 

dependent  on  velocity  source  terms  that  are  linear  (1)  and  non-linear  (2)  with  respect  to  the 
fluctuating  velocity.  This  fluctuating  pressure  can  be  directly  compared  to  the  signals  measured 
during  the  experiments  using  the  wall-mounted  pressure  transducers.  The  integration  in  equation 
(2.10)  can  be  applied  separately  to  the  linear  and  non-linear  terms  in  equation  (2.1 1).  This  allows 
an  estimate  of  the  pressure  fluctuations  associated  with  the  mean-shear,  pms,  and  turbulence- 
turbulence,  ptt,  contributions  to  be  computed  separately  and  compared: 

PmsiV&t)  =  ^  [2AV ifij  4-  B{TJxu2  +  r2tii)  +  2Cr2w2] ,  (2.12a) 

n  m  7 

ptt(r]At)  =  ^  Ai-AyAt[A(utui -ufTiT)  +  -uiu2)  +  C(u2U2 -H2U2)].  (2.12b) 

The  coefficients,  A.  B,  and  C,  arc  defined  in  equation  (2.10),  and  the  mQSE  estimated  velocity 
fluctuations,  Uj,  provide  time-dependent  data  for  the  calculation  of  the  pressure  fluctuations.  Note 
that  the  mean  quantities,  Uj  and  u,u,  ,  are  computed  directly  from  the  PIV  data  and  not  from  the 
mQSE  estimation  results. 

In  (2.12),  the  integration  from  (2.10)  has  been  replaced  by  a  summation  with  weights,  Ax, 
Ay,  and  Az  for  a  trapezoidal  approximation.  In  the  streamwise  and  wall-normal  directions,  the 
spatial  resolution  of  the  PIV  data  dictates  that  Ax  =  Ay  =  0.8  mm.  The  variation  of  the 
quadrupolc  in  the  spanwise  direction  must  be  accounted  for;  however,  velocity  data  was 
measured  (and  therefore  estimated)  in  a  plane  at  only  a  single  spanwise  location  (on  the 
centerline).  Therefore,  for  the  purposes  of  this  calculation,  the  velocity  field  was  assumed  to  be 
homogeneous  in  the  spanwise  direction  such  that  U,(y3)  =  Uj(y’3).  This  is  not  a  valid  assumption 
in  terms  of  the  flow  physics,  but  it  was  used  here  out  of  necessity  to  carry  out  the  calculation  in 
(2.12).  Finally,  the  width  of  the  cavity  was  segmented  into  10  equal  parts  such  that  Az  =  5.08 
mm,  and  the  velocity  field  was  assumed  to  be  the  same  in  each  spanwise  segment. 

Equation  (2.12)  was  used  to  calculate  pms  and  ptt  at  the  floor  locations  where  the  pressure 
fluctuations  were  measured  during  the  experiment.  The  contribution  of  the  linear  and  non-linear 
sources  can  be  summed  to  give  an  estimate  of  the  total  fluctuating  pressure  load,  pk(r|At)  = 
Pms.k(flAt)+p„,k(r|At).  This  allows  a  direct  comparison  between  the  estimated  and  measured 
pressure  loads  for  the  same  instants  in  time. 

2.3.2  Source  Identification  Application 


Figure  9  and  Figure  10  show  the  time  history  of  the  measured  cavity  floor  pressure  loads 
at  Mach  0.19  and  0.73  compared  to  the  estimated  fluctuations  calculated  from  equation  (16).  In 
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each  figure,  time  has  been  normalized  by  the  convection  velocity,  0.57Umf,  and  the  cavity  length, 
L,  such  that  t=  0.57Ujnft/L.  The  results  from  these  two  Mach  numbers  are  representative  of  all 
the  Mach  numbers  studied  here.  At  each  Mach  number,  the  measured  pressure  loads  exhibit  a 
predominant  downstream  propagation  that  occurs  over  a  single  cycle.  This  downstream  motion 
becomes  more  dominant  at  the  higher  Mach  numbers  with  a  more  pronounced  periodicity  and  is 
represented  well  by  the  estimated  pressure  loads.  Compared  to  pms  and  ptt,  only  the  mean-shear 
contribution  to  the  pressure  loads  exhibits  this  same  downstream  propagation  and  periodicity. 
This  observation  suggests  that  the  resonant  features  of  the  cavity  are  primarily  related  to  the 
mean-shear  contribution  to  the  pressure  loads. 


Figure  9:  Time  history  of  the  estimated  wall-pressure  loads  at  Mach  0.19:  (a)  measured 
values  used  to  calculated  the  estimated  velocity  field,  pk(r|At);  (b)  estimated  pressure 
fluctuations,  pk(r|At);  (c)  estimated  contribution  to  the  pressure  from  the  mean-shear  component, 
Pms.k(TlAt);  (d)  estimated  contribution  from  the  non-linear  component,  ptt,k(r]At).  The  time  has 
been  normalized  such  that  t=  0.57Ujnft/L.  Contour  range  is±0.08  psi  and  black  is  negative. 
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Figure  10:  Time  history  of  the  estimated  wall-pressure  loads  at  Mach  0.73:  (a)  measured 
values  used  to  calculated  the  estimated  velocity  field,  pk(hDt);  (b)  estimated  pressure 
fluctuations,  pk(hDt);  (c)  estimated  contribution  to  the  pressure  from  the  mean-shear  component, 
Pms,k(hDt);  (d)  estimated  contribution  from  the  non-linear  component,  ptt-k(hDt).  The  time  has 
been  normalized  such  that  t=  0.57Uinft/L.  Contour  range  is±0.08  psi  and  black  is  negative. 


Figures  7  through  9  show  the  auto-spectral  density  of  the  measured  pressure  loads  on  the 
cavity  floor  compared  to  the  spectra  of  the  estimated  pressure  fluctuations  and  the  contributions 
of  the  mean-shear  and  non-linear  velocity  sources.  The  one-sided  auto-spectral  density  plotted  in 
the  figures,  G((^Af ),  was  calculated  as 

C*(CA/)  =  2(pJ(CA/)pt(CA/))/A/  (2.13) 
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where  Pk(^Af )  is  a  1024  point  FFT  of  the  pressure,  p*  denotes  the  complex  conjugate,  and  <•> 
denotes  ensemble  averaging  over  the  16  blocks  of  data.  Results  for  Mach  0.19,  0.39  and  0.73  arc 
shown  and  represent  the  same  features  found  in  the  results  at  Mach  0.29  and  0.58. 

The  estimates  pi  and  pn  both  under-predict  the  magnitudes  compared  to  the  measured 
pressure  spectra.  This  is  most  likely  related  to  neglecting  the  effect  of  the  end  walls  but  could 
also  be  influenced  by  the  limited  extent  of  the  PIV  data.  Nevertheless,  the  estimated  pressure 
spectra  still  exhibit  the  resonant  features  at  these  measurement  locations. 

The  broadband  levels  (particularly  for  the  non-linear  contribution)  tend  to  be  relatively 
flat  and  arc  often  larger  than  that  of  the  measured  values  at  the  higher  frequencies.  This  is  likely 
a  result  of  the  uncertainty  in  the  PIV  data.  Assuming  that  the  integration  of  equation  (2.12) 
causes  random  error  to  have  little  relevance,  the  bias  uncertainty  in  the  PIV  data  becomes  the 
primary  source  of  error  for  the  pressure  calculation.  Given  a  bias  uncertainty  of  approximately 
0.4  m/s  with  1200  samples  of  the  velocity  in  the  PIV  data,  the  pressure  calculation  could  be 
effected  by  as  much  as  (0.4m/s/1200°5)  •  420Pa  •  s/m  =  4.8Pa  which  is  approximately  107  dB 
referenced  to  20  pPa.  Although  this  is  not  an  exact  treatment  of  the  uncertainty,  it  demonstrates 
that  the  uncertainty  in  the  PIV  data  is  on  the  order  of  the  broadband  levels  in  the  calculated 
pressure  spectra.  This  creates  a  kind  of  background  noise  floor  for  the  estimated  pressure  spectra 
which  would  have  a  larger  effect  at  the  lower  Mach  numbers  when  the  overall  magnitudes  of  the 
pressure  fluctuations  arc  smaller  than  the  higher  Mach  numbers. 

In  some  instances,  the  estimated  pressure  exhibits  a  peak  in  the  spectra  that  is  absent  in 
that  of  the  measured  pressure  or  vice  versa.  There  is  no  clear  explanation  for  this  phenomenon  at 
this  time.  However,  it  is  possibly  the  result  of  the  estimated  pressure  fluctuations  being 
completely  dependent  on  the  velocity  field  without  taking  into  account  the  effects  of  the  cavity 
acoustics,  scattering,  or  other  similar  factors  which  are  known  to  be  important  in  this  type  of 
problem.  Nevertheless,  the  analysis  still  shows  a  clear  distinction  between  the  features  of  the 
linear  and  non-linear  sources  for  the  pressure  loads.  Most  importantly,  the  results  demonstrate 
that  the  mean-shear  or  linear  contribution  to  the  pressure  fluctuations  carries  the  resonant 
features  of  the  measured  pressure  loads.  The  non-linear  contribution  is  primarily  broadband  in 
nature  and  does  not  exhibit  any  of  the  dominant  Rossiter  mode  features.  In  most  eases,  the  mean- 
shear  contribution  even  matches  the  magnitudes  of  the  peaks  in  the  measured  spectra.  Even  at 
Mach  0.19.  where  a  resonance  does  not  clearly  present  itself,  the  mean-shear  contributions 
matches  the  broad  rise  in  energy  at  the  lower  frequencies  below  1 500  Hz.  Therefore,  the  resonant 
features  of  the  cavity  are  reflected  predominately  in  the  mean-shear  contribution  from  Poisson’s 
equation.  These  results  suggest  that  a  successful  cavity  resonance  control  strategy  would  modify 
the  mean-shear  characteristics  of  the  shear-layer.  This  is  in  agreement  with  the  host  of  leading 
edge  devices  used  in  either  a  passive  or  active  sense  to  control  cavity  flows  that  have  been 
studied  by  many  authors  as  reviewed  by  Rowley  &  Williams().  The  fact  that  the  turbulence- 
interaction  features  of  the  shear-layer  are  not  associated  with  the  cavity  resonance  may  help 
determine  why  some  control  strategies  function  more  favorably  than  others. 
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Figure  II:  Auto-spectral  density  of  the  11  floor  mounted  pressure  measurements 
compared  to  the  estimation  results  and  the  contributions  from  the  mean-shear  and  non-linear 
components  for  Mach  0. 1 9.  Pressure  measurement  location,  P|<,  refers  to  the  location  along  the 
cavity  floor.  Magnitude  is  dB/Hz  referenced  to  20  gPa. 
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Figure  12:  Auto-spectral  density  of  the  11  floor  mounted  pressure  measurements 
compared  to  the  estimation  results  and  the  contributions  from  the  mean-shear  and  non-linear 
components  for  Mach  0.39.  Pressure  measurement  location,  Pk,  refers  to  the  location  along  the 
cavity  floor.  Magnitude  is  dB/Hz  referenced  to  20  pPa. 
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Figure  13:  Auto-spectral  density  of  the  1 1  floor  mounted  pressure  measurements 
compared  to  the  estimation  results  and  the  contributions  from  the  mean-shear  and  non-linear 
components  for  Mach  0.73.  Pressure  measurement  location,  P^,  refers  to  the  location  along  the 
cavity  floor.  Magnitude  is  dB/Hz  referenced  to  20  pPa. 
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3  Flow  Physics  of  Successfully  Controlled  Cavity  Flow 

In  the  following  seetion  we  will  eompare  various  aspeets  of  the  velocity  field  for  a 
successfully  controlled  eavity  flow.  The  data  sets  being  used  for  this  application  is  from  a  3- 
dimensionsional  eavity  geometry  sitting  in  a  Maeh  1.5  flow.  The  eontrol  eame  from  an  array  of 
miero-jets  distributed  at  the  leading  edge  of  the  eavity.  These  data  sets  are  both  numerical  and 
experimental.  The  experimental  data  set  was  acquired  by  Prof.  Farrukh  Alvi  and  his  students  at 
Florida  State  University.  While  the  numerical  data  set  was  generated  by  Dr.  Srinivasan 
Arunajatesan  at  Combustion  Research  and  Flow  Technologies  using  their  one  equation  Large 
Eddy  Simulation  model.  In  what  follows  we  will  only  eoneentrate  on  comparing  the  flow 
characteristics  between  the  baseline  (no  eontrol)  and  the  ease  where  the  floe  eontrol  resulted  in 
the  maximum  reduction  of  fluemating  surfaee  pressure.  The  details  of  the  experiments,  numeries 
and  the  effieaey  of  the  eontrol  on  the  surfaee  pressure  ean  be  found  in  Arunajatesan  et  al  (2009), 
Ukeiley  et  al,  (2008)  and  Sheehan  (2007).  It  should  be  noted  that  work  being  presented  in  this 
this  seetion  was  part  of  collaboration  between  this  AFOSR  program  with  other  programs  at 
AFRL/RBAI 

3.1  Cavity  and  Actuation  Description 

The  analysis  being  presented  in  this  seetion  is  for  a  nominal  length  to  depth  ratio  5.6  eavity 
whieh  was  three  dimensional  in  nature  with  slanted  side  walls  and  a  sloping  floor.  A  schematic 
of  the  eavity  is  shown  in  Figure  14. 


Figure  14:  Sehematie  of  the  nonrectangular  eavity  being  studied  in  thepresent  work. 

The  actuation  consisted  eight  400  pm  mierojets  whieh  were  loeated  just  upstream  of  the 
eavity  leading  edge  at  a  distance  of  four  mierojet  diameters  from  the  leading  edge.  The  spacing 
between  the  mierojets  matehes  that  in  the  experimental  work  at  FSU  presented  in  Ukeiley  et  al 
(2008)  and  Sheehan  (2007),  which  was  parametrically  optimized.  A  schematic  of  the  orientation 
and  location  of  the  mierojets  ean  be  seen  in  Figure  15.  The  mierojet  pressures  are  40  psig  and 
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arc  operated  choked.  These  conditions  correspond  to  the  conditions  where  the  microjcts’ 
effectiveness  started  to  saturate. 


Figure  15:  Schematic  of  the  microjet  actuators. 

The  effects  of  the  control  on  the  root  mean  square  surface  pressure  fluctuations  are 
presented  in  Figure  16  for  both  the  experiments  and  the  numerical  simulations.  Clearly  there  is  a 
significant  reduction  in  the  fluctuating  surface  pressure  for  cases.  It  also  should  be  noted  even 
though  there  are  some  differences  between  the  numerical  and  experimental  results  (~  1 0%)  they 
both  show  approximately  the  same  levels  of  reduction.  Additional  it  was  demonstrated  in  the 
Ukeiley  ct  al  (2008)  and  Arunajatesan  ct  al  (2009)  that  there  were  reductions  in  both  the  tonal 
and  broad  band  features  of  the  fluctuating  surface  pressure  spectra. 


Figure  16:  Reductions  of  root  mean  square  surface  pressure  fluctuations  for  both 
experimental  and  numerical  data  sets. 


32 


3.2  Mean  and  Turbulent  Flow  Fields 

Figure  17  displays  contours  of  the  mean  flow  for  both  the  control  and  no  control  cases  for 
both  the  streamwisc  and  wall-normal  components.  From  this  figure  one  can  see  substantial 
changes  in  the  flow  patterns  within  the  cavity.  Examining  the  streamwise  component  one  clearly 
observes  that  the  shear  layer  is  lifted  and  it  does  not  appear  to  have  as  significant  an  impact  with 
the  aft  wall.  The  artifact  of  this  lifting  of  the  shear  layer  can  be  seen  by  examining  the  wall- 
normal  velocity  where  there  is  clearly  less  of  downward  velocity  on  the  aft  wall. 


Mean  WaFNomnal  Velocity 

No  Control 

0 

12  3*56 
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Figure  17:  Mean  flow  velocity  contours  demonstrating  the  effects  of  the  leading  edge 

microjets. 
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Figure  18:  RMS  flow  velocity  contours  demonstrating  the  effects  of  the  leading  edge 

microjets. 

Figure  18  displays  the  contours  of  the  root  mean  square  of  the  fluctuating  velocities  for 
both  the  controlled  and  no  control  scenarios.  Once  again  the  lofting  of  the  shear  layer  is  apparent 
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from  studying  where  the  peak  levels  in  the  streamwise  component  reside.  Also  apparent  from 
these  figures  is  that  the  shear  layer  for  the  controlled  ease  thiekens  quite  a  bit  more  than  the 
uncontrolled  ease  in  the  initial  part  of  its  development.  From  examination  of  the  wall-normal 
rms  contours  there  is  a  reduction  in  the  size  of  the  unsteady  region  at  the  aft  end  of  the  cavity. 
Additionally,  there  is  a  lifting  of  the  most  intense  fluctuations  out  of  the  eavity.  Although  not 
presented  above  the  observations  from  the  numerical  simulations  were  quite  similar  although  the 
effects  on  the  fluctuating  velocity  were  not  as  substantial. 

3.3  Proper  Orthogonal  Decomposition 


An  application  of  the  Proper  Orthogonal  Decomposition  was  earned  out  on  both  the 
experimental  and  numerical  data  sets  discussed  above.  Results  for  the  baseline  eavity  are 
compared  with  the  mierojet  control  eases,  as  these  were  the  only  eases  where  experimental 
flowfield  data  is  available.  The  snapshot  method  (Sirovieh,  1987)  was  used  to  extraet  the  POD 
modes  from  the  flowfield  solution  due  to  its  computational  efficiency  with  highly  spatially 
resolved  data.  The  analysis  discussed  first  only  uses  the  velocity  data  from  the  eavity  midplane. 
The  kinetie  energy  norm  comprising  only  of  u  and  v  (negleeting  w)  is  used  to  eompare  with  the 
modes  obtained  from  PIV  on  a  similar  plane  from  the  experimental  data. 


Figure  19:  Comparsion  of  POD  modes  with  and  without  control. 


The  convergence  of  the  energy  in  the  POD  modes  for  both  the  numerical  and 
experimental  data  sets  for  both  the  baseline  and  controlled  eases  is  presented  in  Figure  19.  An 
interesting  observation  from  these  plots  is  that,  contrary  to  expectations,  the  convergence  of  the 
control  ease  is  better  than  the  baseline  ease.  There  is  more  energy  in  the  first  mode,  whieh  has 
typically  been  interpreted  to  imply  an  increase  in  large-scale  organization  of  the  flow.  As  has 
been  discussed  here  and  elsewhere  (Arunajatesan,  2006),  the  use  of  mierojets  leads  to  a  reduction 
in  the  flow  structure  coherence.  This  behavior  of  the  modes  is  also  observed  experimentally,  as 
shown  in  Figure  19b  where  the  POD  modes  for  the  experimental  data  set  also  show  that  the 
amount  of  energy  in  the  first  mode  is  larger  for  the  controlled  flow  for  the  baseline  even  though 
these  modes  contain  a  higher  percentage  of  the  overall  energy.  However,  one  must  be  careful  in 
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interpreting  these  results,  as  it  will  be  demonstrated  later  that  neglecting  the  vv  in  the  norm  has 
led  to  this  contrary  result,  and  a  full  3-D  POD  may  be  required. 
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Figure  20:  Comparison  of  streamwise  POD  modes  of  baseline  vs.  control  cases. 

A  comparison  of  the  first  three  POD  modes  of  the  streamwise  component  from  both  the 
experiment  and  simulation  results  is  presented  in  Figure  20.  In  general,  the  experimentally  and 
numerically  extracted  POD  modes  show  agreement  with  the  differences  in  spatial  resolution 
being  quite  apparent  though.  For  the  first  mode  the  effect  of  eontrol  is  to  loft  the  shear  layer 
slightly  over  the  eavity,  as  is  shown  qualitatively  by  both  the  simulation  and  experiment. 
However,  the  effect  is  more  pronounced  in  the  experiment.  The  second  and  third  modes  arc 
qualitatively  similar  with  the  third  mode  representing  the  dominant  wavelength  of  structures  that 
is  very  similar  to  that  in  the  experiment.  Further,  the  modes  with  and  without  control  appear  to 
be  more  or  less  similar  in  nature  exeept  in  the  elosc  proximity  to  the  leading  edge  of  the  cavity 
where  the  control  concept  is  applied.  This  inference  reinforces  the  earlier  results  about  the 
turbulence  field  being  only  slightly  different  for  most  of  the  cavity  except  near  the  leading  edge. 

In  an  attempt  to  more  quantitatively  assess  the  similarity  of  the  POD  modes  the  similarity 
function  (Rcmpfcr  ct  al,  1994)  was  calculated  for  the  experimentally  determined  POD  modes  and 
is  presented  in  Figure  21.  This  function  uses  the  fact  that  the  POD  modes  arc  orthonormal  hence 
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the  inner  product  of  two  modes  which  are  the  same  would  be  one.  In  Figure  21  the  inner  product 
of  the  first  20  modes  from  the  baseline  (no  control)  ease  are  taken  with  the  first  20  modes  from 
the  controlled  case.  In  the  plot  (based  on  the  color  code  in  the  legend)  square  which  are  blue 
represent  mode  combinations  where  the  inner  product  is  small  and  the  modes  arc  not  similar.  On 
the  contrary  squares  which  are  red  represent  values  of  the  inner  product  is  approaching  one  and 
hence  these  mode  combinations  represent  modes  which  are  similar  for  both  the  controlled  and 
baseline  cases.  An  example  of  on  such  combination  is  the  fourth  mode  from  the  controlled  flow 
and  the  seventh  mode  from  the  baseline.  The  end  goal  of  studying  this  will  be  to  identify  key 
flow  features  that  exist  in  both  the  baseline  and  controlled  cases  so  more  intelligent  control  can 
be  developed  in  the  future  to  push  the  flow  to  one  dominated  by  these  features. 


Figure  21:  Similarity  function  from  experimental  POD  modes  demonstrating  similarities 
in  the  modes  from  the  controlled  and  baseline  eases. 
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4  Closed-Loop  Active  Control  Studies 


In  this  seetion  we  will  present  results  from  and  application  of  elosed  loop  adaptive 
eontrol  to  a  resonating  eavity  flow.  While  the  application  is  admittedly  at  a  lower  Mach  number 
the  methodologies  developed  here  are  can  be  utilized  at  higher  free  stream  Mach  numbers  if 
actuators  with  enough  authority  can  be  developed.  This  controlled  flow  will  serve  as  base  flow 
in  follow  on  studies  for  us  to  develop  a  better  understanding  of  what  the  eontrol,  which  is  based 
on  reducing  the  surfaee  pressure  fluctuations,  does  to  the  flow. 

In  this  ehapter,  first,  a  fccdbaek  control  methodology  is  developed  for  reducing  the  flow- 
indueed  eavity  oseillation  and  broadband  pressure  fluctuations.  Then,  an  aetive  flow  control 
aetuator  is  designed  and  investigated.  Finally,  adaptive  system  ID  and  eontrol  algorithms  are 
combined  and  implemented  in  real-time  to  achieve  the  flow  control.  In  order  to  aehieve  the 
objectives,  a  MIMO  system  ID  algorithm  are  derived  and  discussed  first.  Then,  the  MIMO 
adaptive  GPC  algorithm  is  described.  This  is  followed  by  a  description  of  the  wind  tunnel 
facilities  and  the  data  processing  methods.  Wind  tunnel  experimental  results  for  both  open-loop 
(baseline)  and  elosed-loop  are  then  presented  and  diseussed.  Finally,  the  eonelusions  and  future 
work  arc  presented. 

4.1  System  Identification  Algorithm 


In  this  seetion,  a  MIMO  system  ID  algorithm  is  developed.  A  linear  system  model  is 
summed  with  the  r  inputs  [ujrxi  and  the  m  outputs  [y]mxi-  For  simplification,  the  order  p  of  the 
feedback  loop  is  assumed  the  same  as  the  order  of  the  forward  path.  At  speeifie  time  index  k,  the 
system  can  be  expressed  as 


y(k)  =  aty(k  - 1)  +  a,  y(k  -  2)  +  •  •  •  +  apy(k  -  p ) 

+finu(k)  +  fru(k  - 1)  +  pzu(k -  2)  +  •  •  •  +  fipu(k  -  p) 

Where 


u(k)  =  [w(£)]rx|  = 

ux(k) 

u2(k) 

y2(k) 

ur(k) 

rxl 

a.  =[«.]  .a,  =\aA 

l  L  1  Jmxm  i  L  *•  Jmxm  p  \_  P  Jmxm 

A  =[A>Lxr>A  '  0 P  ~  [A>]m*r 


Define  the  observer  Markov  parameters 

-  a„  |  A  -  Al.,..,..,,.,,! 


(4.3) 
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and  the  regression  vector 


y(k- •) 


<p{k)  = 


y(k~p) 

"«"(*) 


u(k-p) 


[/n*/?+r*( /?+!)]*  I 


(4.4) 


substituting  Equation  4.3  and  Equation  4.2  into  Equation  4.1  yields  a  matrix  equation  for  the 
filter  outputs 


Furthermore,  the  errors  are  defined  as 


(4.5) 


(4.6) 


Finally,  the  observer  Markov  parameters  4.2  can  be  identified  recursively  by 

0(k  +  \)  =  0(k)-ps(k)q>T(k) 

In  order  to  automatically  update  the  step  size,  choose 

I 

m  = — nr 

o’+HL 


where  a  is  a  small  number  to  avoid  the  singularity  when  is  zero. 


GENERALIZED  PREDICTIVE  CONTROL  ALGORIHTM 


(4.7) 


(4.8) 


The  generalized  predictive  control  (GPC)  algorithm  belongs  to  a  family  of  the  most 
popular  model  predictive  control  (MPC).  The  MPC  algorithm  is  a  feedback  control  method, 
different  choices  of  dynamic  models,  cost  functions  and  constraints  can  generate  different  MPC 
algorithms.  It  was  conceived  near  the  end  of  the  1970s  and  has  been  widely  used  in  industrial 
process  control.  The  methodology  of  MPC  is  represented  in  Figure  22,  where  k  is  the  time  index 

number,  u(k)  are  the  input  sequences,  and  y(k)  are  the  actual  output  sequences.  The  and 
-Vr^'  are  estimated  output  and  reference  signals,  respectively. 

Two  comments  are  made  here  to  describe  all  MPC  algorithms.  First,  at  each  time  step,  a 
specific  cost  function  is  constructed  by  a  series  of  future  control  signals  up  to  u(k+s)  and  a  series 

of  future  error  signals,  which  are  the  differences  between  the  estimated  output  signals  y^  +  j) 
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and  the  reference  signals  >V(*  +  7)  Second,  a  series  of  future  inputs  u(k+j)  are  calculated  by 
minimizing  this  cost  function,  and  only  the  first  input  signal  is  provided  to  the 


Figure  22:  Model  predictive  control  strategy. 


system.  At  the  next  sampling  interval,  new  values  of  the  output  signals  are  obtained,  and  the 
future  control  inputs  are  calculated  again  according  to  the  new  cost  function.  The  same 
computations  are  repeated.  Juang  et  al.  (1997,  2001)  give  the  derivation  of  the  adaptive  MIMO 
GPC  algorithm.  This  algorithm  is  an  effective  control  method  for  systems  with  problems  of  non- 
minimum  phase,  open  loop  unstable  plants  or  lightly  damped  systems.  It  is  also  characterized  by 
good  control  performance  and  high  robustness.  Furthermore,  the  GPC  algorithm  can  deal  with 
the  multi-dimension  case  and  can  easily  be  combined  with  adaptive  algorithms  for  self-tuning 
real-time  applications.  The  problem  of  flow-induced  open  cavity  oscillations  exhibit  several 
theses  issues,  therefore,  the  GPC  is  considered  as  a  potential  candidate  controller. 

Two  modifications  are  made  for  this  algorithm.  First,  a  input  weight  matrix  is  integrated 
into  the  cost  function,  this  control  matrix  can  put  the  penalty  for  each  control  input  signal  and 
further  to  tune  the  performance  for  each  input  channels.  Second,  a  recursive  version  of  GPC  is 
developed  for  real-time  control  application. 
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MIMO  Adaptive  GPC  Model 

In  this  section,  a  MIMO  model  with  HR  structure  is  considered.  A  linear  and  time 
invariant  system  with  r  inputs  [u]rxi  and  m  outputs  [y]mxi,  at  the  time  index  k  can  be  expressed  as 


y(k )  =  or, y(k  - 1)  +  a2y(k  -  2)  +  •  •  •  +  apy{k  -  p) 

+fi0u(k)  +  ftu(k  - 1)  +  fi2u(k  -  2)  +  •  •  •  +  fipu{k  -  p) 

where 


‘«,(^ 

Mj(^) 

»(*) 

v2(A-) 

_ur(k)_ 

rxl 

ar.=[a.]  ,a,=\aA  ,-,a=  \aD\ 

1  L  Jfflxm  ^  l  1  Jmxm  7  7  P  L  P  Jmxm 

A  =[A]m,r>^i  =[A]mx(.’"‘iA  =  [A>]mi<r 


(4.9) 


(4.10) 


Shifting  k  step  ahead  from  the  Equation  (4.9),  the  output  vector  y(k+j)  can  be  derived  as 


y{k  +  j)  =  a^y(k  - 1  )+•••  +  «, V(*  “P  +  0 

+apU).y(*  -/>)+ #,«(*  +  y)  +  ff0(l)u(k  +  j  - 1)+ •  •  • 
+A(y,«(*)+A0V*-1)+-+A<>,iK*-p) 


Where 


a1u_,)a1+a2(>’,) 

L  Jmxm 

[a2,y)] 

a/i'‘naJ+a,(>-n 

Jmxm 

a  ,0,1 

,+a0'1 

P  Jmxm 

1  p-1  p 

r«  w,i 

=  a,°-')an 

L  P  Jmxm 

I  P 

[au’L=(“i"'”a+a''"”) 

[«0'L=(“.u"^+aih’) 
[awL =at‘~"p. 


and  with  initial 


(4.11) 


(4.12) 


40 


(4.13) 


a 


P-\ 


a 


p 


[awL=a 

[AmL=A 

Ia-.i0’L=a,-, 

[^I0'L=a 


/? '*> 

The  quantities  y°  (k=0,l,...)  are  the  impulse  response  sequence  of  the  system. 
Defining  the  following  the  vector  form 

'  ti(k-p)  ^  f  u(k) 

u(k  +  1) 


up(k-p)  = 


yP(k-p)  = 


u(k- p  + 1) 

«(*-!) 

y(k-p)  x 
y(k-p  +  \) 


,uJ+l(k)  = 


/(r.p)xl 


U(k  +  j)){r.U+ 


l))xl 


y{k- 1) 


(4.14) 


/  (m*p)x\ 


the  predictive  index  j=0, 1,2,... ,q,q+l,...,s-l, ,  and 
f  u(k)  }  f  y(k)  N 


«.(*)  = 


u(k  +  \) 

V  u(k  +  s  - 1) 


•y,(k)= 


'  (r*.?)x| 


y(k+ 1) 

y(k  +  s- 1) 


(4.15) 


Finally,  the  predictive  model  for  future  outputs,  ys,  is  obtained,  this  future  outputs 
consists  of  a  weighted  summation  of  future  inputs,  us,  previous  inputs,  up,  and  previous  outputs, 

yP 

ys(k)  =  Tus(k)+  Bup  (k-  p)+  Ayp  (k-  p)  (4.16) 

where 


00  0  -  0 

/r  fio  •••  o 


(4.17) 
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(4.18) 


A 

A-. 

...  fix  - 

A<n 

A-.(1> 

...  /r 

fit" 

n  0-1) 

r  P~\ 

...  pi*-" 

"  an 

■■■  a, 

A  = 

a  (1) 

P 

<v/n 

-  a,,n 

a  u~" 

I  P 

CC'-t" 

- 

1  -l( 

( m*s)x(r*p ) 


( p) 


(4.19) 


MIMO  Adaptive  GPC  Cost  Function 

Assume  the  control  inputs  (present  input  and  future  inputs)  depend  on  the  previous  inputs 
and  output  and  can  be  expressed  as 


=["Lv» 


(/>*{m+r)] 


up{k-p) 

yP(k-p) 


(/?*(m  +  r)]x| 


(4.20) 


Two  potential  cost  functions  are  list  below.  The  first  one  consists  terms  of  future  outputs  and  a 
trace  of  the  feedback  gain  matrix 

J(k)  =  ysT(k)Qys(k)  +  rtr(HTH)  (4.21) 


and  the  second  definition  of  cost  function  based  on  the  total  energy  of  future  outputs  as  well  as 
the  inputs 


J(k)  =  ^(ysT(k)QyXk)  +  uiT(k)RuAk)) 


(4.22) 


The  output  weight  matrix  Q,  input  weight  matrix  R  and  the  control  horizon  s  are 
important  parameters  for  tuning  the  controller.  The  horizon  s  is  usually  chosen  to  be  several 
times  longer  than  the  rise  time  of  the  plant  in  order  to  ensure  a  stable  feedback  controller  (Gibbs 
et  al.  2004).  Also,  if  the  predict  horizon  range  is  from  zero  to  infinity,  the  resulting  controller 
approaches  the  steady-state  linear  quadratic  regulator  (Phan  et  al.  1998). 

MIMO  Adaptive  GPC  Law 

In  order  to  minimize  the  cost  function,  three  approaches  are  considered  as  follows. 

•  Based  on  Equation  (4.21),  the  control  coefficients  can  be  update  using  adaptive 
gradient  algorithm. 

•  Based  on  Equation  (4.22),  the  optimum  solution  can  be  derived.  How'ever,  this 
method  requires  the  calculation  of  a  matrix  inverse,  so  the  computational  complex  is  higher. 

•  Based  on  Equation  (4.22),  the  control  coefficients  can  be  updated  using  an 
adaptive  gradient  algorithm. 
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The  first  approach  is  examined  by  Kegerise  et  al.  (2004).  In  the  next  section,  the  latter 
two  approaches  are  derived. 


MIMO  Adaptive  GPC  Optimum  Solution 


Based  on  the  cost  function  (4.22),  the  goal  is  to  find 


["L 


l(s*r)x[p*(m+r)] 


or 


to 


minimize  the  cost  function.  We  will  show  that  both  minimizing  the  cost  function  (4.22)  with 


respect  to  control  matrix  and  input  vector 

result.  To  simplify  the  expression,  let’s  define 


k(*>], 


(s*r)xl 


will  provide  the  same 


p*(m  +  r  )]xl 


u„(k-p) 

yP(k-P ) 


[p*(m+r)]xl 


(4.24) 


with  some  algebraic  manipulation,  the  optimum  solution  is  obtained 
"<L  =-(T’QT*rYt’Q[B  a][vp] 


(4.24) 


It  is  easy  to  apply  the  optimal  solution  of  the  Equation  (4.24)  on  the  cavity  problem. 
However,  the  matrix  inversion  calculation  has  high  computational  complexity.  Only  if  the 
model  order  is  low  enough,  the  optimal  input  can  be  used  in  real-time  application. 

MIMO  Adaptive  GPC  Recursive  Solution 

To  avoid  calculating  the  inverse  of  the  matrix  in  Equation  (4.24),  the  stochastic  gradient 
descend  method  can  be  used  to  update  the  control  matrix  H  using  the  following  algorithm 

therefore,  the  recursive  solution  is  given  by 

H(k  +  l)=H(k)-p{(TTQT+R]H(k)+TTQ[B  (4.26) 


Since  only  present  r  controls  [u(k)]rxi  are  applied  to  the  system,  only  the  first  r  rows  in 
Equation  (4.26)  are  used 

[h(k  + !,}  =  [h(k»  -  n  {(TrQT*R)[Htk)\  +  TrQ{B  A]}  [vjv,  J  (4  2?) 

first  r  rows 

The  output  weight  matrix  Q,  input  weight  matrix  R  and  the  control  horizon  s  arc  tuning 
for  testing  their  effects  to  the  control  performance. 


4.2  Wind  Tunnel  Experiment  and  Setup 


43 


The  experimental  facilities  and  instruments  used  in  this  study  are  described  in  detail  in 
this  section.  These  devices  consist  of  a  blowdown  wind  tunnel  with  a  test  section  and  cavity 
model,  unsteady  pressure  transducers,  data  acquisition  systems,  and  a  DSP  real-time  control 
system.  Actuator  system  used  in  this  study  is  described  also  in  this  section. 


Wind  Tunnel  Facility 


The  compressible  flow  control  experiments  arc  conducted  in  the  University  of  Florida 
Experimental  Fluid  Dynamics  Laboratory.  A  schematic  of  the  supply  portion  of  the 
compressible  flow  facility  is  shown  in  Figure  23.  This  facility  is  a  pressure-driven  blowdown 
wind  tunnel,  which  allows  for  control  of  the  upstream  stagnation  pressure  but  without 
temperature  control. 

The  compressed  air  is  generated  by  a  Quincy  screw  compressor  (250  psi  maximum 
pressure.  Model  5C447TTDN7039BB).  A  desiccant  dryer  (ZEKS  Model  730HPS90MG)  is  used 
to  remove  the  moisture  and  residual  oil  in  the  compressed  air.  The  flow  conditioning  is 
accomplished  first  by  a  settling  chamber.  The  stagnation  chamber  consists  of  a  254  mm 


Pressure  Valve 


Manual  Valve 


Honeycomb 


Settling  Chamber 


Connect  to 
Test  Section 


Screens 


Figure  23:  Schematic  of  the  wind  tunnel  facility. 
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Figure  24:  Schematic  of  the  test  section  and  the  cavity  model.  Dimensions  are  inches. 


diameter  cast  iron  pipe  supplied  with  the  clean,  dry  compressed  air.  A  computer  controlled 
control  valve  (Fischer  Controls  with  body  type  ET  and  Acuator  Type  667)  is  situated 
approximately  6  meters  upstream  of  the  stagnation  chamber  with  a  76.2  mm  diameter  pipe 
connecting  the  two.  A  flexible  rubber  coupler  is  located  at  the  entrance  of  to  the  stagnation 
chamber  to  minimize  transmitted  vibrations  from  the  supply  line.  The  stagnation  chamber  is 
mounted  on  rubber  vibration  isolation  mounts.  A  honeycomb  and  two  flow  screens  are  located 
at  the  exit  of  the  settling  chamber  and  the  start  of  the  contraction  section,  respectively.  The 
honcyeomb  is  76.2  mm  in  width  (the  cell  is  76.2  mm  long)  with  a  cell  size  of  0.35  mm.  Two 
anti-turbulcnee  screens  spaced  25.4  mm  apart  are  used;  these  screens  have  62%  open  area  and 
use  0. 1  mm  diameter  stainless  steel  wire.  For  the  current  experiment,  the  facility  was  fitted  with 
a  subsonic  nozzle  that  transitions  from  the  254  mm  diameter  circular  cross-section  to  a  50.8  mm 
x  50.8  mm  square  cross-section  linearly  over  a  distance  of  355.6  mm.  The  profile  designed  for 
this  contraction  found  in  previous  work  provides  good  flow  quality  downstream  of  the 
contraction  (Carroll  ct  al.  2004).  The  overall  area  contraction  ratio  from  the  settling  chamber  to 
the  test  section  is  19.6:1.  For  the  present  subsonic  setup,  the  freestream  Mach  number  can  be 
altered  from  approximately  0.1  to  0.7,  and  the  facility  run  times  are  approximately  10  minutes  at 
the  maximum  flow  rate  due  to  the  limited  size  of  the  two  storage  tanks,  each  with  volume  of 
3800  gallons. 

Test  Section  and  Cavity  Model 

A  schematic  of  the  test  section  with  an  integrated  cavity  model  is  shown  in  Figure  24. 
The  origin  of  the  Cartesian  coordinate  system  is  situated  at  the  leading  edge  of  the  cavity  in  the 
mid-plane.  The  test  section  connects  the  subsonic  nozzle  exit  and  the  exhaust  pipe  with  431.8 
mm  long  duet  with  a  50.8  mm  x  50.8  mm  square  cross  section. 
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The  eavity  model  is  contained  inside  this  duct  and  is  a  canonical  rectangular  eavity  with  a 
fixed  length  of  L=1 52.4  mm  and  width  of  W=50.8  mm  and  is  installed  along  the  floor  of  the  test 
seetion.  The  depth  of  the  eavity  model,  D,  can  be  adjusted  continuously  from  0  to  50.8  mm. 
This  meehanism  prov  ides  a  range  of  cavity  length-to-depth  ratios,  L/D,  from  3  to  infinity.  The 
cavity  model  spans  the  width  of  the  test  section  W.  However,  a  small  cavity  width  is  not 
desirable,  beeausc  the  side  wall  boundary  layer  growth  introduces  three-dimensional  effects  in 
the  aft  region  of  the  cavity.  As  a  result,  the  growth  of  the  sidewall  boundary  layers  in  the  test 
seetion  may  result  in  modest  flow  acceleration.  The  boundary  layers  have  not  been  characterized 
in  this  study.  Nevertheless,  the  eavity  geometry  applied  in  this  study  is  consistent  with  previous 
efforts  in  the  literature  (Kegerise  et  al.  2007a, b)  considered  to  be  shallow  and  narrow,  so  two- 
dimensional  longitudinal  modes  will  be  dominant  (Heller  and  Bliss  1975).  Removable,  optical 
quality  plexiglas  windows  with  25.4  mm  thiekness  bound  either  side  of  the  eavity  model  to 
provide  a  full  view  of  the  eavity  and  the  flow  above  it.  The  floor  of  the  cavity  is  also  made  of  14 
mm  thiek  plexiglas  for  optical  access. 

Two  different  wind  tunnel  eavity  eeiling  configurations  arc  available.  The  first  one  is  an 
aluminum  plate  with  25.4  mm  thiekness  that  ean  be  considered  a  rigid-wall  boundary  condition. 
This  boundary  condition  helps  excite  the  eavity  vertical  modes  and  the  “eut-on”  frequencies  of 
the  eavity/duet  configuration  (Rowley  and  Williams  2006).  The  performance  of  this  ceiling  is 
discussed  in  the  next  ehapter. 

In  order  to  simulate  an  unbounded  cavity  flow  encountered  in  practical  bomb-bay 
configurations,  a  flush-mounted  aeoustie  treatment  is  eonstrueted  to  replaee  the  rigid  eeiling 
plate.  The  new  eavity  eeiling  modifies  the  boundary  conditions  of  the  previous  sound  hard 
eeiling.  This  aeoustie  treatment  consists  of  a  porous  metal  laminate  (MKI  BWM  series, 
Dynapore  P/N  408020)  baeked  by  50.8  mm  thiek  bulk  pink  fiberglass  insulation  (Figure  24). 
This  aeoustie  treatment  eovers  the  whole  eavity  mouth  and  extends  1  ineh  upstream  and 
downstream  of  the  leading  edge  and  trailing  edge,  respectively.  This  kind  of  aeoustie  treatment 
reduees  reflections  of  aeoustie  waves.  The  performance  of  this  treatment  is  assessed  in  the  next 
ehapter. 

The  exhaust  flow  is  dumped  to  atmosphere  via  a  5  degree  angle  diffuser  attached  to  the 
rear  of  the  cavity  model  for  pressure  recovery.  A  custom  rcctangular-to-round  transition  pieee  is 
used  to  conneet  the  reetangular  diffuser  to  the  6  ineh  diameter  exhaust  pipe. 

Three  structural  supports  are  used  to  reduee  tunnel  vibrations  (Carroll  et  al.  2004).  Two 
of  these  structural  supports  attach  to  both  sides  of  the  test  seetion  inlet  flange,  and  the  additional 
structural  support  is  installed  to  support  the  iron  exhaust  pipe  (Figure  24). 

Pressure/Temperature  Measurement  Systems 

Stagnation  pressure  and  temperature  are  monitored  during  eaeh  wind  tunnel  run  and 
converted  to  Mach  number  via  the  standard  isentropie  relations  with  an  estimated  uncertainty  of 
±0.01  The  reference  tube  of  the  pressure  transducer  is  eonneeted  to  statie  pressure  port  (shown 
in  Figure  24)  using  0.254  mm  ID  vinyl  tubing  to  measure  the  upstream  statie  pressure  of  the 
eavity.  The  stagnation  and  static  pressures  are  measured  separately  with  Druek  Model  DPI  145 
pressure  transducers  (with  a  quoted  measurement  precision  of  0.05%  of  reading).  The  stagnation 
temperature  is  measured  by  an  OMEGA  thermocouple  (Model  DP80  Series,  with  0.1  deg.  C 
nominal  resolution). 
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Two  pressure  transducers  are  located  in  the  test  section  to  measure  the  pressure 
fluctuations.  The  first  transducer  is  a  flush-mounted  unsteady  Kulite  dynamic  pressure 
transducer  (Model  XT-190-50A)  and  is  an  absolute  transducer  with  a  measured  sensitivity  (2.64 
+/-  0.06)xl0  V/Pa  with  a  nominal  500  kHz  natural  frequency,  3.447x1  O'  Pa  (50  psia)  max 
pressure,  and  is  5  mm  in  diameter.  This  pressure  transducer  is  located  on  the  cavity  floor  (y=-D) 
0.6  inch  upstream  from  the  cavity  real  wall  (x=L),  and  8.89  mm  (z=8.89  mm)  away  from  the 
mid-plane.  This  position  allows  optical  access  from  the  mid-plane  of  cavity  floor  for  flow 
visualization  and  avoids  the  possibility  of  coinciding  with  a  pressure  node  along  the  cavity  floor 
(Rossiter  1964).  The  second  pressure  transducer  is  also  an  Kulite  absolute  transducer  (with 
measured  sensitivity  (5.13  +/-  0.03)xl0  V/Pa  and  nominal  400  kHz  natural  frequency, 
1.724x10  Pa  (25  psia)  max  pressure,  5  mm  in  diameter),  and  it  is  flush  mounted  in  the  tunnel 
side  wall  63.5  mm  downstream  of  the  cavity  as  shown  in  Figure  24.  From  a  series  of  vibration 
impact  tests  performed  in  a  previous  study  (Carroll  ct  al.  2004),  the  results  indicated  that  the 
pressure  transducer  outputs  are  not  affected  by  the  vibration  of  the  structure.  An  experiment  to 
validate  this  hypothesis  is  discussed  in  the  next  Chapter.  Due  to  a  modification  of  the 
experimental  setup,  the  second  pressure  sensor  is  moved  to  the  cavity  floor  (Figure  24)  for  both 
open-loop  control  and  closed-loop  eontrol. 

A  PC  monitors  the  upstream  Mach  number,  stagnation  pressure,  and  stagnation 
temperature,  as  well  as  the  static  pressure.  This  computer  is  also  used  for  remote  pressure  valve 
eontrol  (Figure  23)  in  order  to  control  the  freestream  Mach  number  using  a  PID  controller.  In 
addition,  an  Agilent  E1433A  8-channcl,  16-bit  dynamic  data  acquisition  system  with  built-in 
anti-aliasing  filters  acquires  the  unsteady  pressure  signals  and  communicates  with  the  wind 
tunnel  control  computer  via  TCP/IP  for  synchronization.  The  code  for  both  data  acquisition  and 
remote  pressure  eontrol  output  generation  are  programmed  in  LabVIEW.  The  pressure  sensor 
time-series  data  are  also  collected  for  both  the  baseline  and  controlled  cavity  flows  for  post-test 
analysis. 

Facility  Data  Acquisition  and  Control  Systems 

The  schematic  of  the  eontrol  hardware  setup  is  shown  in  Figure  25.  For  the  real-time 
digital  eontrol  system,  the  voltage  signals  from  the  dynamic  pressure  transducers  are  first  pre¬ 
amplified  and  low-pass  filtered  using  Kemo  Model  VBF  35.  This  filter  has  a  cutoff  range  0.1  Hz 
to  102  kHz,  and  three  filter  shapes  ean  be  used.  Option  41  with  nearly  constant  group  delay 
(linear  phase)  in  the  pass  band  and  40dB/octave  roll-off  rate  is  chosen.  The  cutoff  frequency  is  4 
kHz  for  a  sampling  frequency  of  10.24  kHz.  The  signal  is  then  sampled  with  a  5-channel,  16-bit, 
simultaneous  sampling  ADC  (dSPACE  Model  DS2001). 

The  eontrol  algorithms  are  coded  in  SIMULINK  and  C  code  S-funetions  and  are 
compiled  via  Matlab/Real-Time  Workshop  (RTW).  These  codes  are  uploaded  and  run  on  a 
floating-point  DSP  (dSPACE  DS1006  card  with  AMD  Opteron™  Processor  3.0GHz)  digital 
eontrol  system.  The  DSP  was  also  used  to  collect  input  and  output  data  from  the  DS2001  ADC 
boards  as  well  as  computing  the  control  signal  once  per  time  step.  At  each  iteration,  the 
computed  eontrol  effort  is  converted  to  an  analog  signal  accomplished  using  a  6-ehannel  16-bit 
DAC  (DS2102).  This  signal  is  passed  to  a  reconstruction  filter  (Kemo  Model  VBF  35  with 
identical  settings  to  the  anti-alias  filters)  to  smooth  the  zero-order  hold  signal  from  the  DAC. 
The  output  from  this  filter  is  then  sent  to  a  high-voltage  amplifier  (PCB  Model  790A06)  to 
produce  the  input  signal  for  the  actuator.  The  computer  is  also  able  to  access  the  data  with  the 
dSPACE  system  via  the  Matlab  mlib  software  provided  by  dSPACE  Ine. 
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Figure  25:  Schematic  of  the  control  hardware  setup. 

Actuator  System 

In  order  to  achieve  effective  closed-loop  flow  control,  high  bandwidth  and  powerful 
(high  output)  actuators  are  required.  According  to  Cattafesta  et  al.  (2003),  one  kind  of  actuator  - 
called  “Type  A”  -  has  these  desirable  properties.  Such  actuators  include  piezoelectric  flaps  and 
have  successfully  been  used  for  active  control  of  flow-induced  cavity  oscillations  by  Cattafesta 
(1997)  and  Kegerisc  et  al.  (2002).  Their  results  show  that  the  external  flow  has  no  significant 
influence  of  the  actuator  dynamic  response  over  the  range  of  flow  conditions.  Their  later  work 
(Kegerisc  et  al.  2004;  2007a,b)  also  shows  that  one  bimorph  piezoelectric  flap  actuator  is  capable 
of  suppressing  multiple  discrete  tones  of  the  cavity  flow  if  the  modes  lie  within  the  bandwidth  of 
the  actuator.  Therefore,  the  piezoelectric  bimorph  actuator  is  a  potential  candidate  for  the 
present  cavity  oscillation  problem. 

Another  candidate  actuator  is  the  synthetic  or  zero-net  mass-flux  jet  (Williams  ct  al. 
2000;  Cabell  et  al.  2002;  Rowley  et  al.  2003,  2006;  Caraballo  ct  al.  2003,  2004,  2005;  Dcbiasi  et 
al.  2003,  2004;  Samimy  et  al.  2003,  2004;  Yuan  et  al.  2005).  This  actuator  can  be  used  to  force 
the  flow  via  zero-net-mass  flux  perturbations  through  a  slot  in  the  upstream  wall  of  the  cavity. 
Although  the  actuator  injects  zero-net-mass  through  the  slot  during  one  cycle,  a  non-zero  net 
momentum  flux  is  induced  by  vortices  generated  via  periodic  blowing  and  suction  through  the 
slot. 
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Figure  26:  Bimorph  bender  dise  actuator  in  parallel  operation. 

In  this  research,  a  piezoeleetrie-driven  synthetic  jet  actuator  array  is  designed.  This  type  of 
synthetic  jet  based  actuators  normally  gives  a  larger  bandwidth  than  the  piezoeleetrie  flap  type  of 
actuators.  A  typical  commercial  parallel  operation  bimorph  piezoeleetrie  dise  (APC  Ine.,  PZT5J, 
Part  Number:  P412013T-JB)  is  used  for  this  design.  The  composite  plate  is  a  bimorph 
piezoeleetrie  aetuator,  which  includes  two  piezoeleetrie  patehes  on  upper  and  lower  sides  of  a 
brass  shim  in  parallel  operation  (Figure  26).  The  final  design  of  the  actuator  array  consists  of  5 
single  aetuator  units.  Eaeh  aetuator  unit  contains  one  composite  plate  and  two  reetanglc  orifices 
shown  in  Figure  27.  The  designed  slot  geometries  for  the  aetuator  array  are  shown  in  Figure  28 

Another  advantage  of  this  design  is  that  it  avoids  the  pressure  imbalance  problem  on  the  two 
sides  of  the  diaphragm  during  the  experiment.  Since  the  two  eavitics  on  cither  side  of  a  single 
aetuator  unit  are  vented  to  the  local  statie  pressure,  the  diaphragm  is  not  statically  deflected  when 
the  tunnel  statie  pressure  deviates  from  atmosphere.  The  challenge  is  whether  these  aetuators 
ean  provide  strong  enough  jets  to  alter  the  shear  layer  instabilities  in  a  broad  Maeh  number  range 
and  also  whether  the  aetuators  produee  a  coherent  signal  that  is  sufficient  for  effective  system 
identification  and  control. 
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Figure  27:  Designed  ZNMF  actuator  array.  A)  Operation  plot.  B)  Assembly  diagram  of 
single  unit.  C)  Singe  unit  of  the  actuator.  D)  Actuator  array. 
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Figure  28:  Dimensions  of  the  slot  for  designed  actuator  array. 


Figure  29:  ZNMF  actuator  array  mounted  in  wind  tunnel. 


A  lumped  element  actuator  design  code  (Gallas  et  al.  2003)  was  used  together  with  an 
experimental  trial-and-error  method  to  design  the  single  aetuator  unit.  To  calibrate  this  eompaet 
actuator  array,  the  centerline  jet  velocities  from  each  slot  are  measured  using  eonstant- 
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Figure  30:  Bimorph  3  eenterline  rms  velocities  of  the  single  unit  piezoeleetrie  based 
synthetic  actuator  with  different  excitation  sinusoid  input  signal.  A)  For  side  A.  B)  For  side  B. 

temperature  hotwire  anemometry  (Dantee  CTA  module  90C10  with  straight  general  purpose  1-D 
probe  model  55pl  1  and  straight  short  1-D  probe  support  model  55h20).  A  Parker  3-axis  traverse 
system  is  used  to  position  the  probe  at  the  center  of  actuator  slots.  The  sinusoidal  excitation 
signal  from  the  Agilent  33120A  function  generator  is  fed  to  the  790A06  PCB  power  amplifier 
with  a  constant  gain  of  50  V/V.  The  piezocercmic  discs  arc  driven  at  three  input  voltage  levels: 
50  Vpp,  100  Vpp,  and  150  Vpp,  respectively,  over  a  range  of  sinusoidal  frequencies  from  50  Hz 
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to  2000  Hz  in  steps  of  50  Hz.  Each  bimorph  dise  serves  as  a  wall  between  two  eavities  labeled 
side  A  and  side  B.  The  notation  used  to  identify  eaeh  bimorph  and  its  corresponding  slots  is 
shown  in  Figure  29.  The  rms  velocities  of  the  slots  3A  and  3B  loeated  in  the  eenterline  of  the 
cavity  are  shown  in  Figure  30  as  an  example.  Figure  31  shows  the  simulation  result  calculated 
by  the  LEM  actuator  design  code  and  is  superposed  on  the  experimental  result,  Figure  30.  The 
results  show  that,  the  LEM  aetuator  design  code  provides  a  pretty  aeeurate  rms  velocity 
estimation  of  synthetic  jet  over  a  large  frequency  range  between  50  Hz  and  2000  Hz.  Finally,  the 
measure  input  current  level  to  the  actuator  array  after  the  amplifier  is  measured.  The  results  arc 
shown  in  Figure  32  and  indicate  that  the  input  eurrent  will  saturate  above  136mApp,  whieh 
means  if  the  input  voltage  is  larger  than  100  Vpp,  the  current  to  the  aetuator  will  keep  a  constant 
value.  During  the  closed-loop  experiments,  an  upper  limit  of  150  Vpp  is  used  sinee  the  eurrent 
probe  is  unavailable. 
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Figure  31:  The  comparison  plot  of  the  experiment  and  simulation  result  of  the  actuator 
design  eode  for  bimorph  3.  The  output  is  the  eenterline  rms  velocities  of  the  single  unit 
piezoelectric  based  synthetic  actuator  with  different  excitation  sinusoid  input  signal  for  side  B. 
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Figure  32:  Current  saturation  effects  of  the  amplifier. 
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Figure  33:  Spectrogram  of  pressure  measurement  (ref  20e-6  Pa)  on  the  eavity  floor  with 
aeoustie  treatment  superimposed  with  the  Rossiter  modes  at  L/D=6  (with  a=0.25,  k=0.7). 


Figure  34:  Schematic  of  a  single  periodic  eell  of  the  actuator  jets  and  the  proposed 
interaction  with  the  incoming  boundary  layer. 
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Figure  33  shows  the  spectrogram  of  the  pressure  measurement  on  the  cavity  floor  with 
acoustic  treatment.  The  Rossiter  modes  are  superimposed  on  this  figure.  The  experimental 
details  arc  explained  in  the  next  chapter.  For  this  study,  the  lower  portion  of  the  Mach  number 
range  (from  0.2  to  0.35)  is  our  control  target  as  an  extension  to  previous  work  by  Kegerisc  ct  al. 
(2007a, b).  The  desired  bandwidth  of  the  designed  actuator  should  cover  the  dominant  peaks  of 
Rossiter  mode  2,  3  and  4,  which  is  between  500  Hz  and  1500  Hz.  (Rossiter  mode  1  is  usually 
weaker  compared  to  Rossiter  modes  2,  3  and  4.)  Over  this  frequency  range,  the  designed 
actuator  can  generate  large  disturbances.  In  addition,  the  array  produces  normal  oscillating  jets 
that  seek  to  penetrate  the  boundary  layer,  resulting  in  streamwisc  vortical  structures.  In  essence, 
it  acts  like  a  virtual  vortex  generator.  A  simple  schematic  of  the  actuator  jets  interacting  with  the 
flow  vortical  structures  is  shown  in  Figure  34  The  approach  boundary  layer  contains  spanwisc 
vorticity  in  the  x-y  plane  (the  coordinate  is  shown  in  Figure  24).  By  interacting  with  the  ZNMF 
actuator  jets,  the  2D  shape  of  the  vortical  structures  transform  to  a  3D  shape  with  spanwisc 
vortical  structures.  These  streamwise  vortical  disturbances  seek  to  destroy  the  spanwisc 
coherence  of  the  shear  layer,  and  the  corresponding  Rossiter  modes  are  disrupted  (Arunajatcsan 
et  al.  2003).  Alternatively,  the  introduced  disturbances  may  modify  the  stability  characteristics  of 
the  mean  flow,  so  that  the  main  resonance  peaks  may  not  be  amplified  (Ukeiley  et  al.  2003). 
Unfortunately,  the  flow  interaction  was  not  characterized  in  this  report  and  will  be  addressed  in 
future  work. 

Instead  of  using  one  specific  amplitude  and  one  frequency  in  open-loop  control,  a  closed- 
loop  control  algorithm  is  used  in  this  study  to  exam  the  effects  of  the  disturbance  with  multiple 
amplitudes  and  multiple  frequencies.  Thus,  the  present  actuator  represents  a  hybrid  control 
approach,  in  which  we  seek  to  reduce  both  the  Rossiter  tones  and  the  broadband  spectral  level. 


4.3  Wind  Tunnel  Experimental  Results  and  Discussion 

In  the  present  baseline  (i.c.,  uncontrolled)  experimental  study,  flow-acoustic  resonances 
in  the  test  section  region  and  in  the  cavity  region  arc  examined.  A  schematic  of  the  simplified 
wind  tunnel  model  and  the  cavity  region  of  the  experimental  setup  are  shown  in  Figure  35. 

Data  Analysis  Methods 

The  cavity  and  wind  tunnel  acoustic  modes  can  be  obtained  experimentally  using  two 
approaches.  One  way  is  to  measure  the  output  of  each  unsteady  dynamic  pressure  sensor  for 
different  fixed  freestream  Mach  numbers  and  then  find  the  spectral  peaks  for  each  discrete  Mach 
number.  However,  with  this  method  it  is  difficult  to  track  the  gradual  frequency  changes  with 
Mach  number.  The  other  choice  is  to  record  each  unsteady  pressure  sensor  output  continuously 
as  the  Mach  number  is  increased  gradually  over  the  desired  range.  Then,  a  joint-time  frequency 
analysis  (JTFA)  (Qian  and  Chen  1996)  is  applied  to  these  recorded  pressure  time  series  data. 
JTFA  provides  information  on  the  measurement  in  both  the  time  and  frequency  domains. 
Finally,  the  time  axis  is  converted  to  Mach  number  via  synchronized  measurements  of  the  Mach 
number  versus  time.  Similar  analysis  methods  can  be  found  in  Cattafesta  et  al.  (1998),  Kegerise 
et  al.  (2004),  and  Rowley  et  al.  (2005).  In  this  study,  the  sampling  frequency  for  experimental 
data  collection  is  10.24  kHz  and  the  frequency  resolution  is  5  Hz.  The  cut-off  frequency  of  the 
anti-aliasing  filter  is  4  kHz,  and  500  continuous  blocks  of  time  series  data  are  used  in  the 
analysis.  During  the  experiment,  Mach  number  sweeps  from  0.1  to  0.7  in  about  100  seconds. 
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Figure  35:  Schematic  of  simplified  wind  tunnel  and  cavity  regions  acoustic  resonances 

for  subsonic  flow. 
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Figure  36:  Noise  floor  level  comparison  at  different  discrete  Mach  numbers  with 
acoustic  treatment  at  trailing  edge  floor  of  the  cavity  with  L/D=6. 

Noise  Floor  of  Unsteady  Pressure  Transducers 

The  effective  in-situ  noise  floor  of  the  two  unsteady  pressure  transducers  is  presented  in 
Figure  36.  Each  noise  floor  measurement  is  compared  with  the  spectra  obtained  at  different 
discrete  Mach  numbers  for  the  acoustically  treated  L/D=6  cavity.  Within  the  tested  frequency 
range,  the  signal-to-noise  ratio  (SNR)  is  in  excess  of  30  dB,  which  demonstrates  adequate 
resolution  of  unsteady  pressure  transducers  for  the  present  experiments  despite  their  large  full- 
scale  pressure  ranges. 
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Baseline  Experimental  Results  and  Analysis 

The  unsteady  pressure  transducer  JTFA  measurement  for  the  trailing  edge  floor  of  the 
cavity  is  shown  in  Figure  37.  The  results  illustrate  a  very  elean  flow  field  below  Maeh  0.6. 
Therefore,  the  experimental  Rossiter  modes  Rj  shown  in  JTFA  plot  (Figure  37)  now  follow  the 
estimated  Rossiter  curves.  At  higher  upstream  Maeh  numbers  (M>0.6),  the  experimental 
Rossiter  modes  deviate  slightly  from  the  expected  Rossiter  curves.  This  is  partly  beeause  the 
estimated  eurves  use  the  upstream  static  temperature  to  calculate  the  speed  of  sound.  This 
estimation  docs  not  aceount  for  the  expected  significant  static  temperature  drop  due  to  the  large 
flow  acceleration  near  the  aft  cavity  region  seen  by  Zhuang  et  al.  (2003).  Another  possible 
reason  for  these  deviations  of  the  flow-acoustic  resonanee  eomes  from  the  structural  vibration 
coupling  with  the  Rossiter  modes.  At  high  Mach  numbers  above  0.6,  the  structural  vibrations 
may  cause  a  lock-on  phenomenon  with  the  Rossiter  modes.  For  this  study,  all  experiments  are 
thus  performed  below  M  =  0.6. 

In  conclusion,  the  observed  flow-acoustic  behavior  of  the  acoustically  treated  eavity 
model  behaves  as  expected  below  M  =  0.6  and  is  therefore  suitable  for  application  of  open-loop 
and  closed-loop  flow  control. 
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Figure  37:  Spectrogram  of  pressure  measurement  (ref  20e-6  Pa)  on  the  cavity  floor  with  acoustic 
treatment  superimposed  with  the  Rossiter  modes  at  L/D=6  (with  a=0.25,  k=0.7). 

Open-Loop  Experimental  Results  and  Analysis 

The  open-loop  and  closed-loop  experimental  results  using  the  designed  actuator  array  are 
shown  in  this  section.  Before  the  control  experiments,  measurements  of  the  pressure  sensor  at 
the  surface  of  the  trailing  edge  of  the  cavity  with  the  without  the  aetuator  turned  on  are  shown  in 
Figure  38.  Without  the  upcoming  flow,  the  noise  floor  shows  a  significant  peak  at  660  Hz  and  a 
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Figure  38:  Noise  floor  of  the  unsteady  pressure  level  at  the  surfaee  of  the  trailing  edge  of 
the  cavity  with  and  without  the  actuator  turned  on.  A)  The  exeiting  sinusoidal  input  has 
frequency  1050  Hz  and  amplitude  150  Vpp.  B)  The  exeiting  sinusoidal  input  has  frequency 
1 500  Hz  and  amplitude  1 50  Vpp.  The  peaks  near  600  Hz  and  2000  Hz  are  electronic  noise. 
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Figure  39:  Open-loop  sinusoidal  control  results  for  flow-induced  cavity  oscillations  at 
trailing  edge  floor  of  the  cavity.  A)  At  Mach  number  0.1.  B)  At  Mach  number  0.2.  C)  At  Mach 
number  0.3.  D)  At  Mach  number  0.4.  The  cavity  model  with  6  inch  long  and  L/D=6. 


small  peak  at  2000  Hz.  The  pressure  sensor  can  also  sense  the  acoustic  disturbances  associated 
with  the  excitation  frequency  and  its  harmonics,  and  the  measured  unsteady  pressure  level  can 
reach  115-120  dB.  The  extent  to  which  the  measured  levels  deviate  from  theses  values  with  flow 
on  (considered  below)  indicates  the  relative  impact  of  the  actuator  on  the  unsteady  flow. 


First,  open-loop  active  control  is  explored.  The  purpose  of  the  open-loop  experiments  is 
to  verify  if  the  synthetic  jets  generated  from  the  designed  actuator  array  can  affect  and  control 
the  cavity  flow.  A  parametric  study  for  the  open-loop  control  is  explored  first.  A  sinusoidal 
signal  is  chosen  as  the  excitation  input  with  the  frequency  swept  from  500  Hz  to  1 500  Hz.  The 
open-control  performance  is  best  over  the  frequency  range  1000  Hz  to  1500  Hz,  which 
corresponds  to  the  resonance  frequencies  of  the  actuator  array.  Since  at  the  resonance  frequency, 

the  actuator  array  can  generate  larger  velocity  jet,  and  the  blow  coefficient  ^ c 
increases.  As  a  result,  the  control  effect  increases. 


For  these  open-loop  tests,  the  upstream  flow  Mach  number  is  varied  from  0.1  to  0.4.  For 
illustration  purposes,  results  are  examined  here  for  two  sinusoidal  signals  with  200  Vpp  and 
excitation  frequencies  at  either  1.05  kHz  or  1.5  kHz  to  drive  the  actuator  array.  The  1.05  kHz 
excitation  frequency  is  close  to  the  resonance  frequency  of  the  actuator,  while  the  1.5  kHz 
frequency  lies  between  the  second  and  third  Rossiter  modes.  The  experimental  results  shown  in 
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Figure  39  illustrate  that  this  actuator  array  ean  successfully  reduce  multiple  Rossiter  modes, 
particularly  at  Maeh  number  0.2  and  0.3.  In  addition,  the  pressure  fluctuation  is  mitigated  at  the 
broadband  level  on  the  surface  of  the  cavity  floor  for  all  the  tested  flow  conditions.  However, 
new  peaks  are  generated  by  the  exeitation  frequencies  and  their  harmonies,  especially  at  low 
Mach  number  0.1.  With  increasing  upstream  Maeh  number,  the  unsteady  pressure  level  also 
increases  and  the  effect  of  the  control  is  reduced.  Note  the  synthetic  jets  introduce  temporal  and 
spatial  disturbances  to  modify  the  mean  flow  instabilities  and  destroy  the  eoherenee  structure  in 
spanwise,  respectively.  The  effectiveness  of  the  actuator  seales  with  the  momentum  coefficient, 
whieh  is  inversely  proportional  to  the  square  of  the  freestream  velocity.  So,  as  the  upstream 
Mach  number  increases,  the  synthetic  jets  are  eventually  not  strong  enough  to  penetrate  the 
boundary  layer  and  the  eontrol  effeet  is  reduced.  Future  work  should  perform  detailed 
measurements  to  validate  this  hypothesis. 

The  results  of  the  open-loop  control  suggest  that  this  kind  of  actuator  array  can  generate 
significant  disturbances  not  only  along  the  flow  propagation  direction  but  also  in  the  spanwise 
direction  of  the  cavity.  The  combination  of  these  effects  disrupts  the  Kelvin-Helmholtz  type  of 
convective  instability  waves,  which  are  the  souree  of  the  Rossiter  modes.  As  a  result,  multiple 
resonances  are  redueed  via  active  control.  The  experimental  results  also  show  the  limitation  of 
the  open-loop  eontrol. 

Closed-Loop  Experimental  Results  and  Analysis 

The  open-loop  control  results  suggest  that  this  eompaet  aetuator  array  may  be  effective 
for  adaptive  closed-loop  control.  As  discussed  above,  the  synthetic  jets  add  disturbances  to 
disrupt  the  spanwise  coherence  structure  of  the  shear  layer  and  result  in  a  broadband  reduction  of 
the  oscillations.  However,  at  the  same  time,  the  coherence  between  the  drive  signal  and  the 
unsteady  pressure  transducer  will  be  reduced.  High  coherence  is  considered  essential  for 
accurate  system  identification  methods.  To  exam  the  accuracy  of  the  system  ID  algorithm  with 
the  ehange  of  the  estimated  order,  an  off-line  system  ID  analysis  is  first  performed.  The  nominal 
flow  condition  is  chosen  at  M  =  0.275  (to  match  that  of  Kegerise  et  al.  2007a,b)  with  a  L/D=6 
cavity,  and  two  system  ID  signals,  one  with  white  noise  (broadband  frequency  and  amplitude 
0.29  Vrms  )  and  the  other  with  a  chirp  signal  (amplitude  0.86  Vrms  and  ft=  25  Hz  to  fn  =  2500 
Hz  in  T  =  0.05  see),  are  used  as  a  broadband  excitation  souree  to  identify  the  system.  The 
running  error  variances  the  system  ID  are  shown  in  Figure  40.  It  is  clear  that  the  larger  the 
estimated  order,  p,  the  more  accurate  is  the  system  ID  algorithm.  However,  due  to  the 
limitations  of  the  DSP  hardware,  we  eannot  ehoose  very  large  values  of  the  estimated  order  for 
system  ID  algorithm  on-line. 

One  potential  advantage  of  the  closed-loop  adaptive  control  algorithm  is  that  it  does  not 
rely  exclusively  on  accurate  system  ID.  Figure  41  shows  the  result  of  the  closed-loop  real-time 
adaptive  system  ID  together  with  the  GPC  control  algorithm  for  an  upstream  Mach  number  0.27. 
Based  on  the  above  system  ID  results,  due  to  the  DSP  hardware  limitation,  the  estimated  GPC 
order  and  the  predictive  horizon  are  chosen  as  14  and  6,  respectively.  The  breakdown  voltage  of 
the  actuator  array  restricts  the  exeitation  voltage  level;  therefore,  the  diagonal  element  of  the 
input  weight  penalty  matrix  R  is  chosen  as  0.1.  This  research  represents  an  extension  of 
Kegerise  et  al.  (2007b)  where  the  system  ID  algorithm  and  the  elosed-loop  controller  design 
algorithm  are  used  simultaneously  in  a  real-time  application.  It  is  important  to 
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Figure  40:  Running  error  variance  plot  for  the  system  identification  algorithm.  A)  With 
chirp  signal  as  input.  B)  With  white  noise  signal  as  input.  Upstream  Mach  number  is  0.275, 
L/D=6. 
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Figure  41:  Closed-Loop  active  control  result  for  flow-induced  cavity  oscillations  at 
Mach  0.27  at  the  trailing  edge  floor  of  the  L/D  =6  cavity.  The  control  algorithm  uses  an 
estimated  order  of  14  for  both  the  system  ID  and  GPC  algorithms,  and  the  predictive  horizon  is 
chosen  as  6.  A  chirp  signal  is  used  as  the  system  ID  excitation  source. 
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Figure  42:  Input  signal  of  the  Closed-Loop  active  control  result  for  flow-induced  cavity 
oscillations  at  Mach  0.27  at  the  trailing  edge  floor  of  the  L/D  =6  cavity.  The  control  algorithm 
uses  an  estimated  order  of  14  for  both  the  system  ID  and  GPC  algorithms,  and  the  predictive 
horizon  is  chosen  as  6.  A  chirp  signal  is  used  as  the  system  ID  excitation  source. 

note  that  only  the  system  ID  white  noise  or  chirp  signal  is  used  to  identify  the  open-loop 
dynamics,  and  the  feedback  signal  is  not  used  for  this  purpose.  Clearly,  the  results  show  that  the 
GPC  controller  can  generate  a  series  of  control  signals  to  drive  the  actuator  array  resulting  in 
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significant  reductions  for  the  second,  third,  and  fourth  Rossiter  modes  by  2  dB,  4  dB,  and  5  dB, 
respectively.  In  addition,  the  broadband  background  noise  is  also  reduced  by  this  closed-loop 
controller;  the  OASPL  reduction  is  3  dB.  The  input  signal  is  shown  in  Figure  42. 
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Figure  43:  Sensitivity  function  of  the  closed-loop  control  for  M=0.27  upstream  flow 
condition.  The  estimated  order  is  14,  prediction  horizon  is  6,  and  the  input  weight  R  is  0. 1 . 


The  sensitivity  function  discussed  earlier  is  shown  in  Figure  43.  A  negative  amplitude 
value  indicates  disturbance  attenuation,  while  a  positive  value  indicates  disturbance 
amplification.  The  results  show  that  all  the  points  are  negative,  which  indicates  the  closed-loop 
controller  reduces  the  pressure  fluctuation  power  at  all  frequencies.  The  spillover  phenomenon 
(Rowley  et  al.  2006)  is  not  observed  in  Figure  43.  The  spillover  problem  is  generated  because 
either  the  disturbance  source  and  control  signal  or  the  performance  sensor  output  and  the 
measurement  sensor  output  (feedback  signal)  are  collocated. 

The  Bode’s  integral  formula  is, 

J*  \og\S(ico)\dco  =  Ke(pk )  (4.1) 

k 

where  pk  are  the  unstable  poles  of  the  loop  gain  of  the  elosed-loop  system.  So,  for  a  stable 
system,  any  negative  area  at  the  left  hand  side  of  the  Equation  4.1  must  be  balanced  by  an  equal 
positive  area  at  the  left  hand  side  of  the  Equation  4.1.  However,  for  present  elosed-loop  control 
study,  the  left  hand  side  of  the  Bode’s  integral  formula  is  -38  rad/see,  which  shows  that  Bode’s 
integral  formula  does  not  hold  here.  Since  this  formula  is  valid  for  a  linear  controller,  the 
combination  of  the  adaptive  system  ID  and  controller  is  apparently  nonlinear.  A  more  detailed 
study  is  required  in  the  future  to  validate  this  hypothesis. 

A  parametric  study  of  the  GPC  is  then  studied  by  varying  the  estimated  order  and  the 
predictive  horizon.  Figure  44  and  Figure  45show  that  the  control  effects  improve  with 
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Figure  44:  Unsteady  pressure  level  of  the  closed-loop  control  for  M=0.27,  L/D=6 
upstream  flow  condition  with  varying  estimated  order.  The  prediction  horizon  is  6,  and  the  input 
weight  is  0. 1 .  The  excitation  source  for  the  system  ID  is  a  swept  sine  signal. 
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Figure  45:  Unsteady  pressure  level  of  the  closed-loop  control  for  M=0.27,  L/D=6 
upstream  flow  condition  with  varying  predictive  horizon  s.  The  estimated  order  of  the  system  is 
6,  and  the  input  weight  is  0. 1 .  The  excitation  source  for  the  system  ID  is  a  swept  sine  signal. 
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Figure  46:  Unsteady  pressure  level  comparison  between  the  open-loop  control  and 
closed-loop  control  for  M=0.27  upstream  flow  condition.  The  open-loop  control  uses  a 
sinusoidal  input  signal  at  1150  Hz  forcing  and  150  Vpp  and  the  rms  value  is  53  V.  The  closed- 
loop  control  uses  an  estimated  order  14,  predictive  horizon  6,  input  weight  0.1,  and  the  input  rms 
value  is  43  V. 


increasing  order  and  predictive  horizon. 

The  comparison  between  the  open-loop  and  closed-loop  results  is  shown  in  Figure  46for 
the  same  flow  condition.  Notice  that  the  baseline  measurement  for  a  same  flow  condition  can 
vary  a  little  from  case  to  ease.  The  open-loop  control  uses  a  sinusoidal  input  signal  at  1150  Hz 
forcing  and  150  Vpp  and  the  rms  value  of  the  input  is  53  V.  The  closed-loop  control  uses  the 
estimated  order  14,  the  predictive  horizon  6,  and  the  input  weight  0.1,  and  the  input  rms  value  is 
43  V. 
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5  Development  of  Quantitative  Schlieren  System 


As  mentioned  earlier  in  order  to  be  able  to  improve  the  estimate  of  the  time  dependent 
flow  field  and  calculate  the  sources  of  the  intense  fluctuating  pressures  in  the  cavity  in  a 
compressible  flow  we  will  need  to  determine  quantitative  information  about  the  density  field. 
This  section  details  the  parts  of  this  effort  which  were  devoted  to  developing  a  quantitative 
Schliem  tool.  The  following  section  will  start  with  a  brief  introduction  which  is  followed  by  a 
description  of  the  experimental  setup.  This  is  followed  by  a  sub-section  which  describes  the 
analysis  using  the  Two  Microphone  Method  and  the  quantitative  Schlieren  instrument  to 
investigate  the  acoustic  field  in  a  normal  impedence  tube  where  special  attention  is  provided  to 
the  photodetector,  enabling  the  estimate  of  the  overall  instrument  sensitivity.  Finally  we  will 
presents  the  results  and  discussion  on  experiments  in  the  impedance  tube. 


5.1  Introduction  and  Experimental  Set  Up 

The  schlieren  technique  allows  for  the  visualization  of  large  density  gradients  commonly 
encountered  in  applications  such  as  optical  component  testing  and  high-speed  flows.  The 
technique  proved  indispensable  for  supersonic  flow  imaging  since  World  War  II  when  it  was 
applied  to  weapons  development  (Settles,  2001).  It  has  also  been  used  for  turbulent  flow- 
visualization  that  uses  various  methods  to  produce  density  gradients  such  as  a  thermal  gradient 
(Davis,  1 97 1 )  or  a  gaseous  mixture  e( Al-Ammar  et  al,  1 998). 

Based  on  earlier  work  of  Davis  (1971)  along  with  Wilson  and  Damkevcla  (1970) 
McIntyre  et  al  (1991)  developed  a  modem  schlieren  technique  termed  “optical  deflectometry”. 
The  method  coupled  a  conventional  z-type  schlieren  system  with  photodetectors  possessing  very 
high  frequency  response  to  quantify  density  gradients  in  high-speed  turbulent  shear  flows.  The 
technique  has  been  used  to  quantify  cavity  shear  layers  (Garg  and  Cattafesta,  2001)  and  high¬ 
speed  axisymmetric  jets(Doty  and  McLaughlin,  2005).  Also  of  interest  is  the  effort  by  Garg  and 
Settles  (1998)  to  develop  a  quantitative  focused  schlieren  instrument  which  limits  the  integration 
path  to  approximately  5  mm.  The  tool  developed  here  adapts  the  optical  deflcctometry  technique 
to  measure  the  acoustic  field  in  a  normal-incidence  acoustic  impedance  tube. 

A  schematic  of  the  normal-incidence  acoustic  impedance  tube  is  shown  in  Figure  47.  The  tube 
is  0.724  m  long  with  a  25.4  mm  square  cross  seetion;  it  allows  a  plane  wave  acoustic  field  over  a 
frequency  range  of  0  -  6.7  kHz.  The  walls  of  the  impedance  tube  consist  of  22  mm  thick 
aluminum.  A  22  mm  thick  aluminum  plate  is  placed  at  one  end  of  the  tube  to  simulate  a  sound- 
hard  boundary  condition.  A  compression  driver  (JBL  Pro  2426H)  was  mounted  at  the  opposite 
end  of  the  tube  via  a  circular-to-square  transition  piece.  A  5  kHz  sinusoidal  waveform  is 
generated  using  an  Agilent  33120A  function  generator  and  a  Techron  7540  power  amplifier. 
The  test  section  was  located  at  the  end  of  the  tube  close  to  the  specimen.  In  this  study,  borosilica 
optical  quality  glass  windows  are  located  32.86  mm  from  the  specimen  and  are  170  mm  long, 
25.4  mm  wide,  and  3.175  mm  thick.  The  transition  from  the  test  section  to  the  specimen  is 
smooth  and  care  is  taken  to  minimize  leakage  at  all  joints  in  the  tube. 
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Figure  47:  Schematic  of  normal  incidence  impedance  tube. 


The  setup  of  the  z-type  optical  deflectometcr  is  shown  in  Figure  48Error!  Reference  source 
not  found.,  and  a  list  of  the  labeled  system  components  is  provided  in  the  Appendix.  The  light 
source  consisted  of  a  100  W  tungsten-halogen  lamp  in  air-cooled  aluminum  housing.  An  F2/50.8 
mm  condenser  lens  focused  the  light  on  an  adjustable  rectangular  aperture  set  at  1  mm  by  2  mm, 
which  was  placed  at  the  focal  point  of  the  first  parabolic  mirror  f fl 0/1 00  mm).  The  resulting 
collimated  beam  passes  through  the  test  section,  and  the  second  parabolic  mirror  focuses  the 
beam  on  a  knife  edge  mounted  on  a  1-D  traverse  system  (Newport  Model  ESP100)  for 
nanometer-automated  positioning  during  calibration.  Two  focusing  lenses  were  used  to  project 
the  shlicren  image  on  a  translucent  screen  mounted  on  a  2-D  traverse  (Velmex  Model 
MB4012P40J-S4)  controlled  using  a  Velmex  controller  (Velmex  Model  VXM  1).  A  fiber-optic 
cable  with  a  0.25  mm  diameter  pinhole  aperture  is  mounted  on  the  screen  to  route  the  optical 
signal  to  an  avalanche  photodiode  (Hamamatsu  Model  C5460-01).  A  second  fiber-optic  sensor 
is  used  to  monitor  intensity  variations  of  the  source. 
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Figure  48:  Schematic  of  optical  deflectomcter  setup. 


Two  3.175  mm  diameter  Briiel  and  Kjaer  (Type  4138)  microphones  arc  used  for  simultaneous 
acoustic  pressure  measurements  to  reconstruct  the  standing  plane  wave  pattern  in  the  test  section. 
These  two  microphones,  referred  to  as  Mic  1  and  Mic  2,  arc  placed  t/l=63.8mm  and 

d2  =  44.8  mm  from  the  specimen,  respectively.  They  are  flush  mounted  in  a  rotating  plug  to  the 
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top  of  the  tube,  and  both  microphones  are  connected  to  a  Bruel  and  Kjaer  PULSE  Multi-Analyzer 
System  (Type  3560),  which  served  as  the  power  supply  for  the  microphones.  The  analysis  was 
carried  out  entirely  in  the  frequency  domain  and  ensemble  averaging  was  performed  over  1,000 


periodic  blocks  of  data  with  512  data  points  in  each  block,  thereby  reducing  uncorrelated  noise. 


The  sampling  frequency  of  the  data  acquisition  system  was  set  at  51.2  kHz,  providing  a 
frequency  resolution  of  100  Hz.  The  output  signals  of  the  two  microphones  and  the  two 
photodiodes  were  passed  to  a  16-bit  A/D  converter  (National  Instrument  Model  NI4522)  with 
built-in  anti-aliasing  filters.  The  dc  component  of  the  reference  photodiode  signal  used  to  correct 
for  the  temporal  variations  in  the  light  source  intensity  was  sampled  using  a  HP34970A  data 
acquisition  unit.  The  positioning,  data  acquisition,  and  traverse  systems  were  processed  and 
controlled  using  LabVIEW. 

5.2  Two-microphone  method 

TMM  is  used  to  measure  the  complex  reflection  coefficient  and  the  normal  incident  acoustic 
impedance  (Chung  &  Blaser,  1980  and  Scybcrt  &  Ross,  1977).  This  complex  reflection 
coefficient  completely  specifics  the  acoustic  field  for  planar  modes  and  is  used  to  calculate  the 
density  fluctuation  gradient  along  the  tube  axis.  A  frequency  response  function  Hn  between 
Mic  1  and  Mic  2  is  estimated  by 


(5.1) 


where  Gn  is  the  estimated  cross  spectrum  between  Mic  1  and  Mic  2,  and  Gn  is  the  estimated 

autospectrum  of  Mic  1  using  standard  methods  described  in  reference  0. 

The  complex  reflection  coefficient  is(Scybcrt  &  Ross,  1977) 


(5.2) 


where  j  =  V-T  ,  k  is  the  wave  number  of  the  air  along  the  axis  of  the  tube,  dx  is  the  distance 


from  Mic  1  to  the  specimen,  and  s  =  di-d2  is  the  distance  between  Mic  1  and  Mic  2.  The 

uncertainty  associated  with  this  method  is  described  in  Shultz  et  al  (2007). 

Since  the  TMM  assumes  that  the  sound  field  inside  the  tube  is  planar,  the  time-harmonic 
acoustic  field  is  (Blaekstock,  2000) 


p{d,t)  =  Pm(e'kd  +  Re-'kJ)e^ 
=  P(d)e-Jn* 


(5.3) 


where  Pm  is  the  complex  incident  wave  amplitude.  For  a  linear  iscntropic  process,  pressure  and 
density  fluctuations  are  related  by 


P(d)  =  p(d)c02 


(5.4) 
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where  c0  is  the  isentropic  sound  speed.  Therefore,  the  pressure  gradient  in  thet/  direction  is 

dPjd)  2dp(d) 

dd  0  dd  1 


Substituting  Equation  (5.3)  into  Equation  (5.5)  results  in 

^l  =  PJk(e^-Rejkd)/c02  (5.6) 

The  goal  of  the  approach  in  this  paper  is  to  determine  the  density  fluctuation  gradient  using  the 
non-intrusive  optieal  deflectometer  device. 


5.3  Quantitative  Schlieren  method 


The  static  calibration  establishes  the  relationship  between  the  light  intensity  fluctuation  on  the 
image  screen  and  the  density  gradient  fluctuation  in  the  test  section.  The  angular  defleetion  ed 
is  related  to  the  density  gradient  in  the  d  direction  via  (Settles,  2001) 


kW  dp 
n0  dd 


(5.7) 


where  £  =  2.259x10  4m’/kg  is  the  Gladstone-Dale  constant  for  air  at  standard  conditions,  W  is 


the  width  of  the  test  section  along  the  optical  path,  and  n0  is  the  refractive  index  of  the  medium. 
As  shown  in  Figure  49,  assuming  that  the  knife-edge  is  perpendicular  to  the  direction  of  the 
angular  deflection,  parallel  light  rays  defleeted  by  a  small  angle  sd  result  in  a  displacement  of 
the  light  source  image  by  an  amount 


Ad  — 

where  f2  is  the  focal  length  of  the  second  parabolic  mirror. 


(5.8) 
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Figure  49:  Light  source  image  in  the  plane  of  the  knife  edge  (50%  cut-off). 


Figure  49  shows  the  illuminated  side  of  the  knife-edge,  on  which  a  conjugate  image  of  the 
rectangular  source  appears.  However,  the  knife-edge  cuts  off  a  fraction,  typically  50%,  of  the 
incident  light  to  achieve  the  schlieren  effect.  Assuming  that  the  power  of  the  light  beam  is 
conserved  after  passing  through  the  knife-edge  and  that  diffraction  effects  can  be  neglected,  then 

/  a  -/  a'i)  (5-9) 

where  Ixnm  and  I KE  are  the  intensity  of  the  beam  at  the  screen  and  knife  edge,  respectively,  h 

is  the  height  of  the  source  image  at  the  knife-edge  position,  and  a'  is  the  dimension  of  the  source 
that  is  not  cut  off  by  the  knife  edge.  As  shown  in  Figure  49,  when  parallel  light  rays  are 
refracted  due  to  a  density  gradient  in  the  test  section,  the  light  source  image  shifts  by  an 
amount  Aa .  As  a  result,  the  change  in  the  intensity  on  the  image  screen  is 


A /  = 


Ike^L 


A 


screen 


(1) 


Equations  (5.7)  -  (5.10)  are  combined  to  find  the  optical  sensitivity  of  the  schlieren  instrument, 
defined  as  the  change  in  intensity  per  unit  change  in  the  density  gradient 


A/  kf2W 
op,ical  dp/dd  a\ 


(5.11) 


To  measure  A/,  a  linear  avalanche  photodiode  is  mounted  on  the  image  screen.  This 
photodetector  itself  has  a  sensitivity  Spd=AV/AI  that  converts  intensity  fluctuations  into 

voltage  fluctuations.  The  total  sensitivity  of  the  sensor  is  the  products,^,  =  S  icolS^ .  Clearly,  it 
is  desirable  to  have  a  large  sensitivity,  but  other  considerations,  such  as  sensor  noise  and 
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bandwidth,  should  be  considered  when  designing  the  instrument.  Equation  (5.1 1)  shows  how  to 
increase  Snpljcal ,  namely  by  increasing  f2 .  However,  the  apparent  possibility  of  increasing  IKrem 

and  increasing  the  eutoff  (by  decreasing  a' )  has  important  implications  for  the  choice  of  the 
photodctector.  These  issues  are  considered  below. 

The  device  is  calibrated  to  obtain  the  overall  sensitivity.  The  photodetcctor  is  plaeed  in 
the  1.22x  magnified  screen  image  at  various  locations  in  the  d  direction,  namely  at  thirty 
equally  spaecd  positions  to  cover  one  aeoustie  wavelength  under  study  (5  kHz).  At  eaeh 
location,  the  knife-edge  was  traversed  from  full  to  no  cutoff  to  obtain  the  relationship  between 
the  knife-edge  position,  whieh  is  related  to  the  density  gradient,  andAE.  Note  that  a  statie 
calibration  is  only  required  at  one  point  in  the  test  seetion  for  a  sehlieren  system  with  a  perfectly 
uniform  cut-off.  However,  the  calibration  was  repeated  at  the  30  locations  along  the  test  section 
to  assess  cut-off  non-uniformity  and  improve  the  system  accuracy. 
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Figure  50:  Photodiode  knife-edge  static  calibrations  for  different  locations  in  test  seetion. 


Figure  50  shows  three  typical  static  calibration  curves  for  4th,  15th,  and  24lh  of  the  30  locations  in 
the  test  seetion.  The  zero  knife-edge  position  corresponds  to  a  nominal  50%  cut-off.  The  data 
are  presented  in  normalized  form 


0>  = 


V  -V 

pd  pd,dark 

V  - V 

y  ref  v  ref, dark 


(5.12) 


where  Vpd  is  the  dc  voltage  of  the  photodiode  mounted  on  the  image  screen,  Vref  is  the  dc 


voltage  of  the  reference  photodiode,  and  Vpd  dark  and  Vref  dark  are  their  respective  values  when  no 

light  is  admitted  to  the  sensors.  Equation  (5.12)  thus  accounts  for  any  temporal  variations  in  the 
source  and  also  corrects  for  any  dc  offsets  in  the  photodiodes. 

As  revealed  in  Figure  50,  the  normalized  output  of  the  photodiode  has  a  linear 
relationship  with  the  knife-edge  position  except  near  full  or  no  cutoff  where  diffraction  effects 
are  important.  As  predicted  by  Equation  (5.11),  operating  near  full  eutoff  increases  the 
sensitivity;  however,  the  instrument  is  nonlinear  in  this  regime.  Hence,  the  device  is  operated  in 
the  linear  regime  of  Figure  50. 
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The  effeet  of  non-uniform  eut-off  in  the  image  plane  leads  to  variations  in  the  static 
calibration  curves.  However,  only  the  slopes,  A<t>/Aa,  of  the  statie  calibration  eurves  are 
required  to  calculate  the  density  gradient.  A  calibration  curve  was  found  at  eaeh  of  the  thirty 
locations,  and  the  slope  in  the  linear  region  was  calculated  for  subsequent  data  reduction.  The 
mean  slope  of  the  calibration  curve  was  0.106  mm  with  a  standard  deviation  of  0.005  mm  , 
illustrating  a  variation  in  cutoff  uniformity  of  approximately  5%. 

Equation  (5.12)  defines  O  for  the  undisturbed  ease.  For  the  disturbed  ease  in  the 
presence  of  aeoustie  waves  in  the  tube,  the  time-harmonie  fluctuating  density  gradient  causes  a 
voltage  fluctuation  V '  in  the  photodiode  output,  leading  to 


Ad>  =  d>'-d> 


ref 


r 

-V 

V  ref  .dark 


(5.13) 


Knowledge  of  the  loeal  slope  combined  with  AO  gives  An.  Equations  (5.7)  and  (5.8)  are  then 
used  to  calculate  the  density  gradient  fluctuation.  The  magnitude  and  phase  of  V'  are 
determined  from  the  coherent  power  spectrum  using  the  cross  spectrum  between  the  photodiode 
and  Mie  1. 

The  optical  deflectometer  requires  a  high-sensitivity,  low-noise  photodeteetor.  Three 
common  photodeteetors  include  a  photomultiplier  tube  (PMT),  a  silicon  avalanche  photodiode 
(APD),  or  the  pixels  of  a  CCD  camera."  The  selection  of  the  appropriate  photodeteetor  depends 
on  the  specific  application  and  the  required  signal-to-noise  ratio  of  the  experiment.  Therefore, 
estimating  the  noise  characteristics  of  the  photodeteetor  in  the  optical  configuration  is  important. 
In  addition,  the  noise  model  of  the  photodeteetor  can  be  combined  with  the  quantitative  sehlieren 
device  sensitivity  Slolal  to  calculate  the  minimum  detectable  signal 


=  VnmseISlnla,  (5.14) 

with  dimensions  of  density  gradient.  Unlike  previous  versions  of  the  optical  defleetometer, 
where  large  density  gradients  are  present  (perhaps  due  to  shoek  waves),  the  appropriate  objective 
here  is  to  measure  infinitesimal  density  gradients  corresponding  to  low-amplitude  aeoustie 
waves.  From  Equation  (5.14),  this  is  achieved  both  by  reducing  the  noise  floor  of  the 
photodeteetor  and  maximizing  the  device  sensitivity. 

An  APD  (Hamamatsu  Model  C5460-01)  photodeteetor  was  chosen  for  this  experiment. 
The  APD  possesses  high  bandwidth  and  high  sensitivity.  As  described  in  references  Virgin  et  al 
(2005),  an  APD  becomes  the  best  choice  when  the  background  illumination  level  is  high. 
Conversely,  a  PMT  is  the  proper  choice  for  low  illumination  levels.  In  a  sehlieren  system  that 
uses  a  high  intensity  source,  an  APD  is  more  appropriate.  A  noise  model  for  this  type  of 
photodeteetor  is  derived  below. 

A  linear  photodeteetor  generates  an  internal  current  iin  that  is  proportional  to  the  input 
light  power  L.n 


L  =  SpftL,n  (5.15) 

where  Sp  is  the  anode  radiant  sensitivity,  and  //  is  a  voltage-to-eurrent  conversion  factor.  An 
amplifier  then  boosts  this  input  current  to  generate  an  output  current  ioul  =  G  ■  im .  The  input 
detector  current  consists  of  a  mean  (de)  and  a  fluctuating  (ae)  component 
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(5.16) 


due  to  background  light  power  and  fluctuations  of  the  input  light  power,  respectively. 

A  noise  flow  diagram  is  shown  in  Figure  51.  Two  major  sources  of  noise  are  Johnson 
(thermal)  noise,  /, ,  and  shot  noise,  ishol .  Johnson  noise  is  a  constant  background  noise  and  it  is  a 

function  of  the  Boltzmann’s  constant,  AT,  experimental  ambient  temperature, T ,  the  bandwidth 
of  the  system,  B,  and  the  device  resistance,  R  .  This  quantity  is  expressed  as  an  alternating 
current  and  is  defined  as 


4  KTB 


(5.17) 


Therefore,  Johnson  noise  is  a  constant  background  noise  and  is  independent  of  the  changed  input 
light  power. 

On  the  other  hand,  shot  noise  is  a  current  directly  affected  by  the  incident  light  power, 
and  is  consequently  multiplied  by  the  APD  by  a  gain  factor  G  .  The  total  noise  added  during  this 
multiplication  process  is  from  the  main  noise  sources,  the  background  current,  and  fluctuating 
current,  and  the  dark  current.  In  addition,  a  surface  leakage  current,  which  is  not  affected  by  the 
internal  gain,  is  also  generated  by  the  dark  current.  As  a  result,  shot  noise  can  be  expressed  by 
Virgin  et  al  (2005) 


(5.18) 
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Figure  51:  Measured  noise  sources  and  the  noise  flow. 


where  E  is  the  electron  charge  constant,  F  is  an  excess  noise  factor,  iih  is  surface  leakage 
current,  which  is  the  dark  eurrent  not  subject  to  gain,  and  idg  is  an  internal  current,  which  is  the 
dark  current  subject  to  gain  and  is  described  by 

'd=L+G-'dg  (5.19) 

The  internal  gain  factor,  G  ,  is  a  critical  parameter.  It  is  a  wavelength  dependent  value,  and  the 
exeess  noise  factor  F  is  a  function  of  this  gain  factor.  It  should  be  noted  that  all  the  noise 
sources  generate  white  noise  with  bandwidth  limited  by  the  sampling  frequency  or  the  maximum 
bandwidth  of  the  detector. 

The  final  measured  total  noise  eurrent,  N ,  is  the  Lr  norm  of  Johnson  noise  and  the  shot 

noise, 
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W=V'.J+0! 

I  -  AKTB 

=  ^2 E(i  +  i '+  idg  )BG2F  +  2EidB  +  — 


(5.20) 


The  quantitative  sensitivity  of  the  photodeteetor,  defined  by  the  relationship  between  the 
light  fluctuation  intensity  on  the  image  sereen  and  the  voltage  fluctuation  output  from  the 
photodeteetor,  is  considered.  The  fluctuating  eurrent  and  voltage  after  the  internal  gain  ean  be 
expressed  by 

i'  =  SpvGL'  (5.21) 

and 

~~~SpGL'  (5.22) 

respectively,  where  L'  is  the  input  light  power  fluctuation.  This  quantity  ean  be  related  to  the 
sehlieren  system  by  the  light  intensity  ehange  on  the  image  sereen  as 

L' =  U  ■  Aiaector  (5.23) 

where  i4dcter,or  is  the  area  of  the  deteetor  head  exposed  to  the  light.  Therefore,  the  measured 

voltage  fluctuating  signal  ean  be  rewritten  as 

=  (5.24) 

Thus,  the  quantitative  sensitivity  of  the  photodeteetor,  S£h ,  defined  as  the  rate  of  ehange 

of  measured  photodeteetor  voltage  fluctuating  signal  with  respeet  to  the  input  light  fluctuation 
intensity  ean  be  expressed  as 

s?=Ti =S’GA*«~  <5-25> 


The  signal  to  noise  ratio  (SNR)  for  the  measurement  of  the  light  fluctuation  intensity  is  defined 
by 


SNR  = 


V 

(N/p) 


(5.26) 


In  general,  the  lowest  detectable  limit  is  (F')  =(N / /j)  at  SNR  =  1 .  In  terms  of  the  power, 

this  situation  indicates  the  measured  voltage  fluctuating  signal  has  the  same  power  level  as  the 
noise.  Therefore,  we  arrive  at  a  noise  equivalent  power  (NEP)  expression  from  Equation  (5.22). 


NEP  =  (L')mm 

x  '  min 


g L 


2  E(i  +  i '+  iJs)BGlF  +  2Ei<kB  + 


AKTB 

R 


mspg 


(5.27) 
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Total  quantitative  sensitivity  for  both  the  schlieren  system  and  the  photodctector  device, 
is  defined  as  the  ratio  between  the  measured  voltage  fluctuation  from  the  deteetor  and  the  density 
fluctuation  gradient 


V' 

r*  total  T 

q  _7^Y 


dx 


(5.28) 


whieh  ean  be  separated  into  two  parts 


gtotal  _ ^5c/j  S  ^  _ 


M 

Vdp) 

dx J 


\r_ 

A/ 


(5.29) 


Substituting  Equations  (5.1 1)  and  (5.25)  for  S^h  and  S* ,  respectively,  the  full  expression  of 
the  total  quantitative  sensitivity  of  the  optieal  dcflectometer  is  expressed  by 

tot  a I  r» Sch  np 

<7  _  ^<7 

.  ano 

JM-.v 

r  5<Tf<*n.max 

ano 

where  Kcreen.m „  =  s PGI screen  max  Ad'Uc,or  is  the  maximum  voltage  the  deteetor  ean  measure  during 

the  experiment.  Now,  the  minimum  detectable  signal  (MDS)  of  the  density  fluctuation  gradient 
ean  be  ealeulated  from  Equation  (5.28) 


screen  max 


[SpGAd«ec,or] 


(5.30) 


(F’)  .  =S“ 

V  /mtn  <7 


total 


dp 
dx  j 


(5.31) 


Similar  to  the  previous  NEP  operation,  the  power  generated  from  the  MDS  of  the 
deteetor  has  the  same  power  level  of  the  total  noise  generated  by  the  detector  itself.  Thus,  the 
MDS  is  obtained  by  substituting  Equaiton  (5.20)  into  Equation  (5.31) 


MDS  = 


V 


2  E(i  +  i '+ )BG2F  +  2Ei(kB  + 


4  KTB 
R 


total 


(5.32) 


The  final  MDS  ean  be  expressed  by  substituting  Equation  (5.30)  into  Equation  (5.32) 
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AKTR 

J2£(i  +  i  ’+  idg  )BG2F  +  2EidB  + 


4ATfl 


4A.TA 


(5.33) 


Schiieren  System 


Equation  (5.33)  can  be  further  simplified  by  noting  that  /'  is  much  less  than  /  .and  /,  and  that  i 
is  normally  much  larger  than  and  ids .  For  the  present  experiment,  the  shot  noise  due  to  the 
background  light  level  dominates  over  the  Johnson  noise,  resulting  in 


(5.34) 


kWf,V 


screen,  max 


5.4  Application  in  Plane  Wave  Impedance  Tube 

Two  experimental  methods  are  employed  simultaneously  to  investigate  the  density 
fluctuation  gradient  in  the  tube:  the  non-intrusivc  optical  defleetometer  and  the  intrusive  TMM  to 
validate.  The  optical  defleetometer  technique  can  capture  one-dimensional  acoustic  fields  over  a 
range  of  aeoustie  intensities.  Several  eases  of  varied  sound  pressure  level  (SPL)  are 
implemented  with  the  range  from  94  dB  to  128  dB,  as  shown  in  Figure  52.  The  distance,  plotted 
on  the  x-axis,  and  the  density  fluctuation  gradient,  plotted  on  the  y-axis,  are  normalized  by  the 
aeoustie  wavelength  and  Pmjk/c02,  respectively.  For  each  data  point  shown,  a  number  of 

voltage  readings  are  taken  from  the  photodeteetor  and  averaged  in  the  frequency  domain  to 
reduee  the  random  background  noise  of  the  detector.  In  the  high  sound  pressure  level  eases 
(>100  dB),  only  1000  blocks  of  readings  are  taken,  but  as  the  SPL  decreases  below  100  dB,  5000 
bloeks  of  data  are  taken  to  reduee  the  random  error  of  the  measurement. 
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(e)  Input  SPL  is  101  dB. 
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(b)  Input  SPL  is  113dB. 
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(d)  Input  SPL  is  94dB. 


Figure  52:  Normalized  Density  gradient  of  1-D  pressure  wave  vs.  non-dimensional 
distance  as  found  using  optical  deflectometer  and  TMM. 


At  lOOdB  SPL,  the  density  fluctuation  gradients  from  the  two  methods  are  very  similar.  In  the 
ease  of  the  lower  SPL  (<100  dB),  the  results  obtained  with  the  optical  deflectometer  show  a 
similar  trend  with  the  TMM.  The  random  errors  in  the  density  gradient  field  are  large  due  to  the 
low  eoherenee  between  the  microphone  and  the  photodiode.  The  error  bars  are  mostly  found  to 
fall  within  each  other’s  range.  At  128dB  SPL,  the  two  methods  achieved  their  lowest 
uncertainties.  This  ease  is  selected  as  a  testbed  for  density  gradient  comparison. 

The  result  of  the  MDS  in  Equation  (5.33)  depends  on  the  accuracy  of  the  noise  model, 
N ,  and  the  total  sensitivity,  S'"'al .  The  sensitivity  is  compared  with  the  optieal  deflectometer 

and  TMM  results  as  shown  in  Figure  53.  The  three  eurves  mateh  well,  indicating  the  total 
sensitivity  model  provides  a  good  estimate  of  the  density  fluctuation  gradient  in  the  test  section. 
The  minimum  detectable  density  fluctuation  gradient  of  the  APD  in  the  d  direction  of  the  tube 
is  estimated  and  shown  in  Figure  54  according  to  Equation  (5.33).  The  maximum  light  input 

power  is  estimated  as  6  /jW  at  the  fully  open  knife-edge  position.  Notice  that  the  quantity  /  in 
Equation  (5.33)  is  a  constant  anode  current  generated  due  to  the  input  light  power.  In  a  lower 
range  of  input  light  power  (<10  ),  thermal  noise  dominates.  Since  the  shot  noise  is  directly 

proportional  to  the  input  light  power,  with  increasing  light  levels,  shot  noise  eventually  beeomes 
the  dominant  noise  term. 
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Figure  53:  Density  gradient  calculated  from  different  method  for  the  SPL=128  dB  case. 


Figure  54:  General  minimum  detectable  signal. 


The  experimental  power  spectrum  of  the  photodetector  for  the  128  dB  ease  is  shown  in  Figure 
55.  Two  peaks  are  clearly  seen  in  the  d  direction  at  5  kHz  due  to  the  speaker’s  single  tone 
excitation  signal.  Other  than  the  peaks  at  5  kHz  and  de,  the  graph  shows  that  the  variation  of  the 
noise  power  in  the  d  direction  is  very  low,  implying  a  relatively  uniform  light  intensity  along 
the  image  screen.  It  should  be  noted  that  the  experiment  used  the  half  knife-edge  cut-off 
location,  but  to  obtain  higher  sensitivity,  a  different  cut-off  could  be  used  (Settles,  2001) 
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Figure  55:  Power  spectrum  of  the  photodiode. 
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Figure  56:  Noise  equivalent  power. 


The  analytical  and  estimated  noise  equivalent  power  arc  shown  in  Figure  56  The  plot  shows  that 
the  noise  model  gives  a  good  estimate  of  the  noise  floor  (within  25%)..  The  error  between  the 
two  curves  comes  from  the  estimated  errors  of  fluctuating  current,  i',  the  dark  current  subject  to 
a  gain,  id.  ,  and  surface  leakage  current,  ids . 


The  MDS  is  estimated  using  the  noise  model,  Equation  (5.20),  and  it  is  superimposed  on  the 
density  fluctuation  gradient  plot  (Figure  53).  The  signal  to  noise  ratio,  which  is  defined  by  the 
density  fluctuation  gradient  measured  by  the  TMM  and  the  noise  level,  is  shown  in  Error! 
Reference  source  not  found..  It  is  clear  that  the  SNR  is  very  large  for  this  case. 
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Figure  57:  Signal  to  noise  ratio  of  the  TMM  and  the  MDS. 
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6  Summary  and  Future  Work 


In  this  study  we  have  put  together  many  pieces  that  can  be  used  for  adaptive  control  of 
open  cavities  and  understand  its  effects  on  the  flow  physies.  These  include;  the  development  of  a 
methodology  to  estimate  the  time  dependent  velocity  field  (mQSE),  the  application  of  Poisson’s 
equation  to  investigate  the  sources  of  the  fluctuating  surface  pressure,  investigation  of  the  effects 
of  the  turbulent  flow  structure  of  a  successfully  controlled  flow,  the  development  and  application 
of  an  adaptive  flow  control  methodology  for  the  for  reducing  the  fluctuating  surface  pressure  and 
the  development  of  a  quantitative  Sehlieren  technique.  While  at  first  glance  this  may  seem  like  a 
disjoint  list  it  will  provide  the  tools  and  fundamental  understanding  so  that  in  a  future  effort  we 
will  be  able  to  successfully  control  supersonic  flow  over  an  open  cavity  with  a  detailed 
understanding  of  it  effects  on  the  flow  physies.  In  what  follows  we  will  present  a  summary  of 
the  key  findings  from  each  of  the  aspects  of  this  project. 


The  mQSE  provides  a  method  for  estimating  the  temporal  dynamics  of  a  flow  field  using 
a  low-dimensional  model  of  the  flow.  For  cavity  flows,  the  mQSE  allows  the  large  scale 
dynamics  of  the  flow  to  be  predicted  from  surface  pressure  fluctuations  in  the  cavity.  This  work 
has  demonstrated  the  application  of  the  mQSE  as  an  analysis  tool  to  study  subsonic  resonating 
cavity  flow.  Application  of  the  mQSE  to  the  velocity  showed  that  the  flow  at  Mach  0. 19  and  0.29 
was  dominated  by  a  circulation  region  occupying  the  cavity  that  oscillates  back  and  forth  in  time. 
Any  formation  of  or  downstream  propagation  of  vortex  structures  was  intermittent.  At  Mach 
0.39,  0.58  and  0.73,  the  formation  and  downstream  propagation  of  vortices  was  clearly  observed 
and  was  always  aeeompanied  by  a  negative  pressure  fluctuation  on  the  cavity  floor.  The 
dominance  of  a  particular  Rossiter  mode  was  also  clear  based  on  the  number  of  vortex  structures 
in  the  shear  layer.  The  aft-wall  pressure  loads  were  shown  to  be  at  a  minimum  when  a  vortex 
structure  was  either  passing  over  or  impinging  on  the  aft  wall.  The  pressure  would  then  increase 
to  a  maximum  as  the  shear  layer  was  pulled  into  the  cavity  by  the  circulation  of  the  next 
approaching  vortex.  Although  not  presented  above,  the  investigators  have  also  demonstrated  that 
the  mQSE  can  be  applied  to  Sehlieren  images  allowing  the  temporal  dynamics  of  the  density 
gradients  to  be  examined  in  relationship  to  the  velocity  field  with  respect  to  the  wall-pressure 
loads.  This  lays  the  ground  work  for  a  follow  on  study  where  the  investigators  will  use  the 
quantitative  Sehlieren  developed  here  to  investigate  both  the  velocity  and  density  fields 
synchronously. 

The  pressure  fluctuations  on  the  floor  of  an  L/D  =  6.0  rectangular  cavity  were  estimated 
by  integrating  Poisson’s  equation  with  flow  data.  The  mQSE  was  utilized  to  estimate  the 
temporal  dependence  of  the  velocity  field  so  that  a  time-resolved  estimate  of  the  pressure  loads 
could  be  calculated.  By  decomposing  the  velocity  into  mean  and  fluctuation  components, 
Poisson’s  equation  shows  that  the  pressure  depends  on  sources  that  are  linear  (mean-shear)  and 
non-linear  (turbulence-turbulence)  with  respect  to  the  fluctuating  velocity.  This  allowed  an 
estimate  of  the  surface  pressure  loads  related  to  these  two  sources  to  be  determined  separately 
and  compared  in  light  of  the  experimentally  measured  pressure  fluctuations.  A  comparison  of 
the  estimated  pressure  loads  with  the  measured  values  revealed  that  the  estimation  provided  a 
reasonable  approximation  to  the  actual  pressure  loads.  It  was  suggested  that  the  minor 
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discrepancies  that  were  noticed  could  be  the  result  of  ignoring  various  acoustic  effects  like 
scattering  and  interference  while  assuming  the  flow  was  completely  incompressible  (a  necessary 
assumption  for  Poisson’s  equation).  Using  the  estimation  procedure,  the  contributions  of  the 
mean-shear  and  non-linear  sources  were  able  to  be  examined  separated.  The  non-linear 
contribution  was  found  to  be  broadband  in  nature  without  any  clear  tonal  features.  However,  the 
mean-shear  contribution  matched  both  the  frequencies  and  magnitudes  of  the  spectral  peaks 
associated  with  the  Rossiter  modes  that  were  observed  in  the  measured  pressure  signals. 
Therefore,  only  the  mean-shear  contribution  to  the  pressure  loads  was  found  to  be  associated 
with  the  resonant  features  of  the  cavity  flow.  The  next  step  will  be  reformulate  the  pressure 
integration  procedure  to  use  a  compressible  form  of  Poisson’s  equation  and  the  estimated 
velocity  and  density  information  from  the  mQSE  procedure. 

A  novel  piezoclectric-driven  synthetic  jet  actuator  array  is  designed  for  this  research. 
The  resulting  actuator  produces  high  velocities  (above  70  m/s)  at  the  center  of  the  orifice  as  well 
as  a  large  bandwidth  (from  500  Hz  to  1500  Hz)  which  is  sufficient  to  control  the  Rossiter  modes 
of  interest  at  low  subsonic  Mach  numbers.  This  actuator  array  produces  normal  zero-net  mass- 
flux  jets  that  seek  to  penetrate  the  boundary  layer,  resulting  in  streamwise  vortical  structures. 
These  streamwise  vortical  disturbances  destroy  the  spanwise  coherence  of  the  shear  layer,  and 
the  corresponding  Rossiter  modes  are  disrupted.  Alternatively,  the  introduced  disturbances 
modify  the  stability  characteristics  of  the  mean  flow,  so  that  the  main  resonance  peaks  may  not 
be  amplified.  Next,  a  MIMO  system  ID  IlR-based  algorithm  is  developed  based  on  the  structure 
inferred  from  the  global  model.  This  system  ID  algorithm  combined  with  a  GPC  algorithm  is 
applied  to  a  validation  vibration  beam  problem  to  demonstrate  its  capabilities.  The  control 
achieves  ~20  dB  reduction  at  the  single  resonance  peak  and  -9  dB  reduction  of  the  integrated 
vibration  levels.  Finally,  this  control  methodology  is  extended  and  applied  to  subsonic  cavity 
oscillations  for  on-line  adaptive  identification  and  control.  Open-loop  active  control  uses  a 
sinusoidal  signal  with  200  Vpp  and  an  excitation  frequency  of  either  1.05  kHz  or  1.5  kHz,  which 
are  detuned  from  the  Rossiter  frequencies,  to  drive  the  actuator  array.  Multiple  Rossiter  modes 
and  the  broadband  level  at  the  surface  of  the  trailing  edge  floor  are  reduced.  However,  when  the 
upstream  Mach  number  increases  (greater  than  Mach  number  0.4),  the  effects  of  the  synthetic 
jets  from  this  actuator  are  gradually  reduced.  Adaptive  closed-loop  control  is  then  applied  for  an 
upstream  Mach  number  of  0.27;  the  estimated  GPC  order  is  14  and  the  predictive  horizon  is  6. 
To  avoid  saturation  in  the  control  signal,  the  input  weight  penalty  is  chosen  as  0.1.  The  GPC 
controller  can  generate  a  series  of  control  signals  to  drive  the  actuator  array  resulting  in  dB 
reduction  for  the  second,  third,  and  fouth  Rossiter  modes  by  2  dB.  4  dB,  and  5  dB,  respectively. 
In  addition,  the  broadband  background  noise  is  also  reduced  by  this  closed-loop  controller  (i.c., 
the  OASPL  reduction  is  3  dB).  However,  unlike  previously  reported  closed-loop  cavity  results,  a 
spillover  phenomenon  is  not  observed  in  the  closed-loop  control  result.  As  discussed  in  Chapter 
1,  the  spillover  problem  is  generated  by  a  linear  controller  because  the  disturbance  source  and 
control  signal  or  the  performance  sensor  output  and  the  measurement  sensor  output  (feedback 
signal)  are  collocated.  The  nonlinear  nature  of  the  adaptive  system  may  be  responsible  for  this 
effect. 


An  non-intrusive  optical  deflectometry  technique  was  shown  to  be  effective  for 
visualizing  and  measuring  the  acoustic  field  in  a  plane  wave  tube.  The  density  fluctuation 
gradient  obtained  with  the  optical  deflectometer  was  compared  to  that  obtained  with  the  standard 
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two-microphone  method.  The  comparison  indicated  that  the  instrument  can  successfully  detect 
density  fluctuation  gradients  above  a  certain  sound  pressure  level,  successfully  constructing  the 
acoustic  field  for  sound  pressure  levels  ranging  from  101  dB  to  128dB,  suggesting  the  device 
could  be  applied  to  study  the  acoustic  field  in  jet  engine  nacelle  liners.  The  sensitivities  of  the 
sehlieren  setup  and  the  photodeteetor  were  derived  separately.  The  total  quantitative  sensitivity 
was  then  obtained  by  combining  the  sehlieren  and  photodetector  sensitivities.  The  total 
sensitivity  model  was  used  to  compare  to  the  two-mierophone  method,  and  the  results  show  a 
good  match.  This  sensitivity  model  can  be  used  as  a  cost  function  to  optimize  optical  systems 
for  various  applications.  The  noise  model,  noise  equivalent  power,  and  minimum  detectable 
signal  arc  defined  and  derived.  The  simulation  result  of  the  NEP  calculated  by  the  analytical 
noise  model  is  compared  to  the  experimental  measurement.  The  two  methods  provided  similar 
results  of  the  photodetector  noise  floor.  Finally,  the  noise  floor  of  the  device  was  calculated,  and 
the  signal  to  noise  ratio  was  estimated  using  the  TMM.  With  the  sensitivity  of  the  current  system 
the  next  step  will  be  apply  this  in  a  supersonic  cavity  flow  field  with  the  eventual  goal  of 
sampling  simultaneously  with  the  surface  pressure  so  that  it  can  be  used  in  the  mQSE  to  estimate 
the  time  dependent  properties. 
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